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The purpose of this study was to investigate the interaction of 
plane elastic waves with a thin, hollow, cylindrical shell embedded in an 
elastic medium. 

The cylindrical shell is considered to he elastic, isotropic, 
homogeneous, and of infinite length. It is surrounded hy an elastic, iso- 
tropic, and homogeneous medium whose motions conform to the ordinary theory 
of elasticity. A plane stress wave, either dilatational or shear with a 
step variation in time, whose wave front travels in a direction perpendic- 
ular to the cylinder axis, envelops the shell. Later, a Duhamel integral 
is used to study other wave shapes for the incident stress. 

The response of the shell is studied hy expressing the two com- 
ponents of displacement, radial and tangential, in terms of Fourier series, 
each term of which is called a mode. The equations of motion of the shell 
in vacuo are derived from expressions giving the strain and kinetic ener- 
gies due to generalized external forces . Forces on the shell result from 
the stresses in the medium at the boundary. Stresses in the medium are 
taken to he the sum of the stresses due to the incoming stress wave expressed 
in terms of Fourier series whose coefficients are known, and those due to 
the reflected and diffracted effects expressed in terms of a pair of dis- 
placement potentials representing waves diverging from the axis of the 
shell. 

The equations to he solved consist of two pairs of coupled integro- 
differential equations in the generalized coordinates of the shell and the 



2 



displacement p.tentials. By use of a digital computer they are solved mode- 
wise by a step-by-step iterative integration technique known as the Newmark 
Beta Method, with which values for the potential functions^ and the accelera- 
tions^ velocities^ and displacements of the shell are determined. Stresses 
in the shell are found from the displacements^ and the values of the potential 
functions permit determination of stresses for any point in the medium. 

Although the equations are written to include an infinite number 
of modes, only the first three modes are considered in detail. The computed 
solution is compared with values obtained from a series expansion of the 
equations, which is valid for short times, and with the static solution based 
on the theory of elasticity to which the general solution should approach 
asymptotically; In addition, the results of two particular problems are 
compared with results given in another study. 

NiJimerical solutions are obtained to determine the effect of the 
several parameters which describe the relative physical properties of the 
shell and medium. Results are presented in tabular and graphical form. 
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CHAPTER I 
INTRODUCTION 



1.1 General Remarks 

The problem of designing underground protective structures to 
resist the effects of nuclear weapons has become increasingly important in 
recent years with the development of modern weapons whose destructive capacity 
is overwhelming. Engineers in this field are hampered to a great extent by a 
lack of theoretical information on how structures in media such as soil or 
rock behave when subjected to dynamic loads. Even for static loads alone^ 
much of the design practice today is of a semi-empirical nature. 

When a nuclear explosion occurs^ stress waves are transmitted 
through the air and ground. How are they transmitted and how are they modified 
by the presence of a structure embedded in the medium? How does the structure 
respond? 

The purpose of this report is to study one aspect of the problem^ 
the interaction between the medium and structure. 

1.2 Statement of Problem 

The problem considered here consists of analyzing the elastic 
response of a hollow cylindrical shell (tunnel lining) embedded in an elastic 
medium when subjected to an incident plane stress wave traveling in a direction 
perpendicular to the axis of the shell. Some questions with which this problem 
may be associated are: Do tunnel linings in contact with rock afford a measure 

of protection significantly higher than an unlined opening? What magnitude 
and time variation of displacement^ velocity^ and acceleration would equipment 
mounted within such a structure be subjected to? How are stress waves within 
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the medium modified in the vicinity of the shell? This study was conducted 
in an attempt to find some qualitative and quantitative answers to these 
questions within the limitations imposed by the assumptions noted below. 

lo5 .Basic Assumptions 

The cylindrical shell is considered to be of infinite lengthy and 
is embedded in an elastic medium of infinite extent in all directions. A 
plane stress wave whose front travels in a direction perpendicular to the 
cylinder axis envelops the shell. Strains parallel to the axis in both the 
medium and shell are assumed to vanish; thus, since each cross section of the 
shell is exactly similar to every other^ the problem is reduced to one of 
plane strain. 

In the mathematical development of the problem certain basic 
assumptions were made^ the most important of which are given here with a few 
explanatory remarks: 

(1) The medium is considered to be homogeneous ^ isotropic^ and 
linearly elastic. This implies that the ordinary theory of stress wave 
propagation applies. In view of the non -homogeneous^ non-isotropic^ and 
non-elastic characteristics of most materials encountered in nature^ this 
is a severe limitation; however^ current theories of stress propagation 
through such media have not advanced to the stage where this limitation can 
be readily overcome. In the case of some rocks^ though^ this assumption 
may be justified. 

(2) The material in the shell is also considered to be homogeneous ^ 
isotropic^ and linearly elastic. Generally speaking^ this assumption is valid 
for values of stress below the so-called proportional limit of materials 
commonly used. In addition ^ the thickness of the shell relative to its radius 
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is assumed small; this permits expression of all stress components of the 
shell in terms of a function which describes the deflection of its middle 
surface o This deflection must satisfy a linear partial differential equation 
with the appropriate boundary conditions c 

( 3 ) The incident stress wave considered is either a plane dilata- 
tional or a plane distortional (shear) wavoo Under actual conditions^ both 
waves are propagated with the shear lagging the dilatational wave- The 
combined effect for elastic conditions may be determined through the principle 
of superposition c 

(4) The radial and tangential particle velocities of the medium 
at the boundary are equal to that of the shell » This is the continuity 
relation insuring that the shell and the medium are in contact with no 
relative slip occurring at the boundary c 

( 5 ) Any additional mass within the shell is assumed to be dis- 
tributed symmetrically about the axis- The significance of this assumption 
is found in the development of the equations of motion to account for any 
additional mass located within the shell. 

Other assumptions a re presented in the formal development of the 
mathematical expressions used to describe the behavior of shell and medium. 

1 o 4 Method of Approach 

The two components of shell displacement^ radial and tangential, 
are written in terms of Fourier series from which expressions giving the 
strain ener^ and kinetic energy of the shell in vacuo are derived. The 
equations of motion are derived from Lagrange's equations in terms of the 
displacement functions and forces acting on the shell. Forces on the shell 
result from the stresses in the medium at the boundary. Stresses in the 



medium are taken to "be the sum of the stresses due to the incoming stress wave 
expressed in terms of Fourier series whose coefficients are known^ and those 
due to the reflected and diffracted effects expressed in terms of a pair of 
displacement potentials representing waves diverging from the axis of the 
shell. The form of these potentials as derived hy Lamb (4) ^ is in terms of 
sine and cosine series^ the nature and treatment of which has been studied 
by Paul ( 6 )p 

The equations of motion^ described as a pair of coupled integro- 
differential equations^ are solved modewise using a numerical technique 
known as the Newmark Beta Method (5) which permits determination of the 
coefficients of the potentials and values of acceleration^ velocity^ and 
displacement of the shell. 

The solution obtained is compared with values obtained from a series 
expansion of the equations^ which is valid for short times ^ and with the static 
solution based on the theory of elasticity which the machine solution should 
approach asymptotically . 

1.5 Previous Work 

The problem stated above has been the subject of a recent report 
by Baron (l)^ whose analysis consists of first solving for the displacements 
caused at the boundary of an unlined cylindrical cavity subjected to a plane 
stress wave. This is done through an integral transform approach,, the solu- 
tion of the transformed equations being expressed using Hankel functions. 

The evaluation of the inverse transform is accomplished only with great 
computational effort. Values of displacements obtained are then used as 
influence coefficients in determining the displacements of the shell. 



Numbers in parenthesis refer to the corresponding entry in the Bibliography. 
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Solutions obtained for two different sets of parameters by Baron are compared 
in Section 6 with results obtained by the method of solution outlined in 
this report o 



stress wave incident on an unlined cylindrical cavity in an elastic medium. 
The reflected and diffracted waves are described in terms of displacement 
potentials which represent outgoing shear and dilatational waves. A method 
was developed for determining values of these potentials^ and a similar 
method is used in this report. 

lo6 Notation 

Notation is defined throughout the text where it first appears; 
however, the following list summarizes the main uses of certain symbols. In 
discussions of special topics other meanings may be ascribed to the symbols, 
at which time they will be redefined. 



The study by Paul (6) consisted of analyzing the effect of a plane 



A 



Area of cross section of the shell, per unit length. 




Coefficients of Fourier series for stresses in the 
medium. 




Coefficients of Fourier series for total stresses 
in the medium. 




Weighting factors. 




Generalized coordinates for the displacements of 
the medium at the boundary. 



a 



ns 



ns 



Generalized coordinates for displacements of the 
shell o 




Velocities of wave propagation, dilatational and 
shear, respectively «. 



d 



Distance from neutral axis of shell to its outermost 
fiber. 



6 



E, E 


= Moduli of elasticity for the medium and shelly 
respectively 0 Bar over the symbol refers to the 
plane strain modulus. 


e 


= Volumetric strain. 


F, G 


= Generalized coordinates of the displacement 
potentials . 


I 


= Moment of inertia of the shelly per unit length. 


k 


= Parameter which relates to the shape of a time 
dependent stress wave. 


k 

r 


= P/r^ ratio of radii. 


k 

c 


= ratio of velocities. 


m 


= Mass of the shell per unit surface area. 


m* 


= Additional mass within the shell. 


n 


= Mode number. 


Q ^ Q 


= Generalized forces. 


R 


= Radius of the shell. 


\ ’ Sn 


= Multiplying factors. 


r, 0 


= Polar coordinates. 


T 


= Kinetic energy of the shell. 


T' 


= Kinetic energy of additional mass within the shell. 


t 


= Thickness of the shell j also, time. 


u 


= Strain energy of the shell. 


u 


= Displacement vector. 


U, V 


= Components of the displacement vector in the radial 
and tangential directions^ respectively. 


U V 

s 


= Radial and tangential displacement components of 
the shell. 


u ^ u 
x" y 


= Components of the displacement vector in the x and y 
directions^ respectively . 


z 


= Rectangular coordinates. 
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a 

n 


a E 

sn s 

~ ' 'r "g — ’ non-dimensionalized radial component of 

p displacement. 


Pn 


b E 
s n s 

^ ^ y non-dimens ionalized tangential component 

p of displacement. 




= Strains along the x and y axes, respectively. 


^9 


= Circumferential strain. 


C 


= A variable of integration. 




E 

= — , a parameter. 


\ 


= ^ , a parameter. 

K 


\ 


= , a parameter. 


\ 


= ^ ^ a parameter. 




= ^ , a parameter. 


'H 

V 


= ^ pstameter. 


9 ^ r 


= Polar coordinates. 




= Position angle of the incident wave. 


k 


= Curvature of the shell. 


[1 


= Lame^ constants. 


V 


= Poisson* s ratio. 


V 


V 

1-v 




= Arguments of the functions F and G, respectively. 


P 


= Mass density of the medium. 


Ps 


= Mass density of the shell. 


a , a 
xx"' yy 


= Normal components of stress parallel to the 
X and y axis, respectively. 


a 

xy 


= Shearing stress component. 



f i 

4-mi 
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°9Q 


= Radial and tangential normal stresses in polar 
coordinates 0 Additional subscripts s and m^ 
when used^ refer to the shell and medium. 


°r9 


= Shearing stress in polar coordinates . 


^sb 


= Bending stress in the shell. 


0 

smax 


= Maximum stress in the shell. 


0 

P 


= Amplitude of incident dilatational wave. 


0 

s 


= Amplitude of incident shear wave. 


T 


= Non-dimensionalized time in the case of the 
incident dilatational wave. 




= Non-dimensionalized time in the case of the 
incident shear wave. 


9, ^1^ 


= Displacement potentials^ dilatational and shear ^ 
respectively . 
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CHAPTER II 
BASIC EQUATIONS 



2.1 Equations for the Medium 

The differential equation of motion of a particle in a homogeneous, 
isotropic, and linearly-elastic medium in terms of its displacement vector u 
is given by Kolsky (j) in the form 

------ 

(\ + 2p) V V • u - pV X (v X u) = p — 5- (2-1) 

dr 

where X and p are the Lame constants, and p is the density of the medium. The 
vector u may be expressed as the sum of two displacements, the of 

a scalar potential and the curl of a vector potential 

u = Vcp + ^xf (2-2) 



Here cp is a potential giving rise to an irrotational displacement and 'f' a 
potential leading to an equivoluminal displacement. 

If the wave equations (2-3) are satisfied and u expressed as in 
(2-2), the equations of motion (2-l) are automatically satisfied. 



2 

V 9 



1 

c^ dt^ 



V 'a 



2 



dt^ 



(2-3) 



where the velocities of wave propagation are 




(2-^) 



for the dilatational and shear waves ^ respectively » 
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Since the prohlem is essentially a two '-dimensional one wherein only the 
rotation about the cylinder axis is considered^, the expression for the 
propagation of shear waves can be written in terms of a scalar potential 
function 



t ^ 1-1 . vith (2-5a) 

Cg S'b 

cp and t 3iTe functions of x^ y, and to The components of the displacement 
vector u can be expressed as 



u = 

X 



u 






(2-5) 



Stress components are 



0 


11 

CD 


+ 


2[J.€ 


XX 








a 


= ur 






xy 


xy 




a 


(D 

11 


+ 


2il€ 


yy 









(2-6) 



where e represents volumetric strain_, and strain along the x and y axis 

respectively^ and T the shearing strain » The strain components for small 

xy 

displacement are 



X 



X 



Su 

^ 

y ^ 



(2-6a) 



xy 



5u 5u 

+ _JL 
by bx 



In terms of the potential functions^ the equations for the stress components 



become 



11 



XX 



mJ ^ J|) * 



= (J- (2 



xy 



a 



:^2 x2 , .2 , 

d cp a If a t \ 

' 3x^ ay^ 



yy 



mJi ^ J|) ^ ^.(0 - gi> 



(2-7) 



In polar coordinates 

2 _ ^ 1 ^ ^ 

^ ■ ar= " " r= 38® 

and the wave equations expressed in polar coordinates are 



. 2 ^ r ar 2 .„2 2 .2 

or r oO ot 



1 at ^ 

ar" " 38^ ' at" 



(2-8) 



Displacements in polar coordinates become 



u = 



^ 1 ^ 

^r r bo 



V = 



1 Sep ^ 

r bo br 



(2-9) 



where u and v represent radial and tangential components of displacement, 
respectively o 



Stress components may be written as 



a 

rr 



(\ + 2[i) 



~a\ 1 ^ _ i_ at' 


\ 


^ 1 a^9 1 at a^t ' 


^ 2 r brbO ” 2 ^ 

-dr r 


+ •” 
r 


^ r ^g2 r ^ ardS_ 



a 



rO 




2 

r aro0 



2 _ ^ ^ 1 ^ 

2 ae . 2 2 . _2 

r dr r as 



1 ^ 
r dr_ 



9Q 



— ^ ^+2u ^ 



dep 1 a^ 1 at a^t 

dr r ^q2 r ^ drdS 



+ \ 




1 a^t . A- M 

r drdQ 2 bO 



( 2 - 10 ) 
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where ^^P^^sent respectively the radial^ shear and circum- 

ferential (hoop) stresses in the mediumc 

In the following development of the equations^ the incoming wave^ 
considered a step function of time^ is either a dilatational or shear wave. 

The step wave was chosen as a convenient form from which^ through the principle 
of superposition using Duhamel's integral^ the effects of any time dependent 
wave may be approximated. It is shown by Kolsky ( 5 ) that for either the 
incoming dilatational or shear wave^ the reflected wave from a plane body must 
have both dilatational and shear components in order to satisfy boundary con- 
ditions. The same holds in the present case« 

2 . 11 Dilatational Wave 

The orientation of the incoming dilatational (compressional) wave 
with respect to the spatial coordinates and time are shown in Figo 1. Owing 
to symmetry about the x axis^ the radial stress^ hoop stress^ and the radial 
component of displacement are even functions of 0^ while the shear stress and 
tangential component of displacement are odd functions of G. 

In the expressions previously given ^ the dilatational potential may 
now be separated into two parts 

cp = cp. + cp . 

^ ^in ^out 

where 

cp. = - a ^ (x+c,t) 

^in p 1 / 

represents the potential of the incoming step wave^ and cp represents the 

ouu 

potential of the outgoing dilatational wave. 

The displacement and velocity of a particle behind the plane wave 
front can be derived from the potential of the incoming wave^ with notation 



defined in Fig. 1 
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u = - a "Y-' ix-i-c t) 
X p '1 



u 



a 

pc. 



The velocities in polar coordinates become 



u = “ — cos ^ 

pc^ 



V = —2— sin Q 

pc. 



( 2 - 11 ) 



( 2 - 12 ) 



The radial and tangential displacements at the boundary r = R due to the 
incoming plane stress wave are expanded in Fourier series as a function of 
time as 



u(R^t,e) = ^ a^(R,t) cos n0 
n.zO 

( 2 - 15 ) 

00 

v(R^t^0) = ^ b^(R;t) sin nG 
n = l 



and the velocities as 

00 

u(R;t;0) = ^ a^(R^t) cos nG 
n=0 

00 ( 2 - 1 ^) 

v(R^t^Q) = ^ 'b^(R^t) sin nG 
n=l 

The Fourier coefficients are determined as a function of G^^ a measure of the 
degree of envelopment of the shell hy the incoming stress wave (Fig* l)^ by 



evaluation of the integrals 



cos 0 cos nO d0 
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2o 



a = 
n 



b = 
n 



«pci Jq 



2a 

-J£- 

JtpC 



/„ 



sin 6 sin nO d.0 



1 '^0 



The indicated integrations are performed using the orthogonality properties 
of the sine and cosine functions to get the following result 



n(in) JTpc 






sin 

sin 20 



sin(n-l)9 sin(n+l)^ 
- + 



n-1 



n+1 



sin 20. 



n 



(in) npc-i 



sin(n-l)0^ sin(n+l)0^ 



n-1 



n+1 



n = 0 
n = 1 



n = 2, 5, o . . 

( 2 - 15 ) 



n = 1 

n = 2^ 5^ • • • 



The subscript (in) has been added to designate the coefficients due to the 
incident wave. 0^ is a time dependent variable which varies from 0 to it 
during the transit of the incoming wave across the cavity^ after which it 
becomes a constant equal to jc. 

The stresses in the medium behind the wave front are as shown in 
— * 

Figo 1^ where ^ ratio relating the stress in a direction perpendicular to 
the direction of wave propagation^ is derived from the assumption that there 
is no strain behind the wave front in the direction parallel to the wave front. 
It is related to Poisson* s ratio^ of the medium in the following way 



* The bar over the symbol does not indicate a vector quantity. 
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V = rT^ 

The stresses behind the wave front are in polar coordinates 

a = - 0 (cos^Q + V sin^S) 
rr p 

“re = °p 



(2-16) 



= -a (sin 0 + v cos 0) 

UU D 



( 2 - 17 ) 



Stresses in the medium as a function of radius and time are now 
expressed in terms of Fourier series 



rr 



^ A^(r,t) cos nS 



n=0 



re 



(r,t,0) = ^ B^(r,t) sin n© 



(2-18) 



n=l 



00 

O00(i’,t,0) = ^ C^(r,t) cos ne 



n=0 

and the coefficients determined in the same manner as used previously for 
velocities are 

_ .1 _TT. 

n = 0 



(l+v)e^ + (~^) sin 29^ 



A \ = “ 
n(in) 2n 



+ (l+v) sin 20^ + sin n = 2 

(l-v)sin(n-2)0 (l-v)sin(n+2)0 2(l-fv) sin nO 



( 2 - 19 ) 



n-2 



n+2 



n = 1,5, 5, 



B /. 



ap(l-v) 



'n(in) 23T 



sin 40 

sin (n-2 )0^ sin(n+2)0^ 



n = 2 



n-2 



n+2 



n = 1.3A.5, 



n = 0 



C 

n 



(in) 



a 

2jt 



{ — — 

(1+v)0;L ■ 

< (v-l)s^ + (l+v) sin 20 sin 40^ 



n = 2 



(v“l)sin(n-2)0^ (v-l)sin(n+2)0^ 2(l+v)sin nG^ 

' ..... -{- — — . — + ' ■ ' ' 

n-2 n+2 n 

n = 1,5, 4, 5,. 

\ 

where the subscript (in) again refers to the incident stress wavec 

As noted earlier^ the reflected and diffracted wave is represented 
as a diverging or outgoing cylindrical wave from a line source whose origin is 
the axis of the cylinder^ and includes both dilatational and shear components. 
The outgoing potential functions are expressed as infinite series ^ each term 
of which will henceforth be called a mode 
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out 



f (r.t) cos nO 
n ' ^ / 



n=0 






out 



^ sin n0 

n=i 



( 2 - 20 ) 



It is shown in the appendix that the coefficients are of the form 



f 

n 



n 

^1 




00 




” cosh U-, ) cosh nu-, du, 
c, 1 11 



•— cosh Ug) cosh nUg du^ 



( 2 - 21 ) 



where F*^(t - — cosh u_ ) and G*^(t - — 
= 1 1 =2 
derivative of the respective functions, 



cosh u^) represent the nth order 
n being equal to the mode number 



considered. 
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To express the velocity components of the medium at the boundary 
due to the outgoing waves as in equations (2-14-), one differentiates 
equations (2-9) with respect to time and substitutes the expressions given 
above for the potential functions. This results in the following for the 
velocity coefficients 



a 

n 



(out) 



00 

^ / F^^^(|-, ) cosh u, cosh nu, du^ 

n+1 J '1' 1 11 



00 

■~ Y - J ^2 ^^^2 ^^2 



^n(out) n+1 



n 



(- 1 ) 

n+1 



n 




( 2 - 22 ) 



where and represent the arguments of the respective functions 
The subscript (out) refers to the outgoing wave. 

The stresses in the medium due to the outgoing waves are also 
expressed in series form as in equations (2-l8); the coefficients being 
found by substitution into equations (2-10) 



(- 1 )^ 

^n(out) ~ n+2 



/ ) cosh nu + 2 cosh^u-, ] du 

\Jr\ i i M- i i 



00 



(2-23) 
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B 



n(out) 



00 

= ^ ^ ) sinh 2u, sinh nu-, du^ 

n+2 1 1 n 



00 

- ^ cosh 2u^ cosh nu^ du^ 
n+2 ^2'^ 2 2 2 



"2 ° 



-ilili. 



"nCout) ^n+2 



oo 

f F"*^a ) 

'dp, 



2 

cosh nu^ [— - 2 sinh u, ] du^ 
1 |j, I-" 1 



00 

L ^“2 ■'“2 '^“e 



<=2 



2.12 Shear Wave 

The orientation of the incoming plane shear wave with respect to 
the spatial coordinates and time are shown in Fig. 2. In this case the shear 
stress and tangential component of displacement must be even functions of 9 
while the radial stress^ hoop stress^ and radial component of displacement 
are odd functions of 9. 

The incoming wave expressed in terms of a shear potential function 



IS 



= a "X- (x + c^t) 
in s s 2 

from which the velocity of a particle behind the wave front in polar 
coordinates may be derived as 

o 

. s 

u = sin Q 

PCo 



(2-24) 



(2-25) 



V = 



pc. 



cos 9 



The stresses in the medium behind the wave front are as shown in Figo 2. In 
terms of polar coordinates they are 
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a = a sin 20 
rr s 



0^=0 cos 20 
r0 s 

<^69 = -^s 



(2-26) 



Velocities and stresses due to the incoming plane stress wave as functions 

of radial distance and time are expanded in terms of Fourier series in a 

manner similar to that of the dilatational wave^ except that because of the 

difference in symmetry^ the sine and cosine terms are interchanged » The 

coefficients are also found in a similar manner and are given hereo The 

coefficients for the velocity series are 

sin 20^ 



n 



(in) Jtpc, 



X 



0 



1 2 

sin(n-l)0 sin(n+l)0 



n-1 



n+1 



sin 0 



1 

sin 20, 



n(in) jTpc^ 



X < 0, + 



sin(n-l)0^ sin(n+l)0^ 



n-1 



n+1 



n = 1 



n = 2,3;.. 



n = 0 



n = 1 



n = 2,3; 
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and the coefficients for the stress series are 

" sin he 

c h --T+ 

A / . ^ ~ X i 

n(in) Tt I . / ^\r^ ./ 

^ ' sin(n-2)0^ sin(n+2)9^ 

n-2 n+2 



n = 2 

n = 1, 5,^.5, • • ■ 



B 

n 



(in) 



sin 20, 



sin 40, 



X S 0, + 



sin(n-2)0^ sin(n+2)0^ 



n-2 



n+2 



n = 0 



n = 2 

(2-28) 

n = 



20 



*^n(in) ■^n(in) 



0^ is a time dependent variable similar to that for the dilatational wave 
except that now it is a function of the position of the incoming shear wave, 
which travels with velocity Cg. 

The outgoing potential functions are expressed modewise as 



9 , 



out 



= f^(r,t) sin n0 



n=l 






out 



n=0 
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where the coefficients are of the form 



f 

n 



00 




“ — cosh ) 

C-, 1 



- — cosh u^) 

=2 2 



cosh nu^ du^ 



cosh nu^ du^ 



(2-30) 



The terms have the same significance as in the case of the dilatational wave. 
Modal coefficients are derived as "before. The coefficients for the velocity 
terms are 



n(out ) 



00 

= - ^ J' F^’*’^(|^) cosh u^ cosh nu^ du^ 



00 



( 2 - 31 ) 
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^n(out) 



n+1 



^n+2 



(I ) sinh u sinh nu, du. 



+ 



LDl 

n+1 




ii^) cosh Ug 



cosh nUg dUg 



and the coefficients for the stress terms are 



i ( out ) 



00 

= f cosh nu^tj + 2 cosh^ Uj^] dUj^ 

Cl 0 

00 



C2 -0 



B 



n(out ) 



00 

= “ ) sinh 2 Ut sinh nu, du, 

n+2 J_ '1^ 1 1 i 



- cosh nug du. 



Cg ~ 0 



(2-32) 



00 



c 



n(out) 



n+2 

Cn 



^n+2(. ) nu-, [— - 2 sinh^ ] du^ 

1 1 jj, 1 1 



(-l)’^U 

n+2 



00 




sinh 2 u 2 



sinh nu^ du^ 



2.2 Equations for the Shell 

In the following derivations of the equations of motion for the 
shelly it is assumed that the thickness of the shell is small in relation to 
its radius. This assumption permits description of the behavior of the shell 
in terms of its middle surface^ and effectively concentrates all of the mass 
in a line thickness. This should be kept in mind when interpreting the results = 
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The effects of shear and rotatory inertia are not considered, since they 

affect only the very high modes, which are not considered in this study. 

The deformation of any point (d) on the shell can then be described 

by specifying the components of displacement (Fig. ^) : 

= radial displacement, positive outwards 

V = tangential displacement, positive in the direction 
of increasing 9. 

The strain energy of the shell per unit length can be expressed as 



U = 



E I 
s 

2 



2Jt g ^ 2Jt 

f (k)^Rd0 + -^ f 

J Q '•J 0 



(2-55) 



where Es 
I 
A 
k 

R 



"plane strain" mo(iulus of elasticity of the shell 
moment of inertia per unit length of shell 
area of the cross section per unit length of shell 
change in curvature of the cross section 
extensional strain in the 9 direction 
radius of the cylinder 



Flugge (2) gives the expressions for the curvature and extensional strain 
terms as 



n2 

T o u 

1 / S \ 

= -p (— ^ + U ) 



and 



1 , 

‘e • E (55- * “5) 



correct to first order terms. The strain energy equation in terms of the 
conponents of displacement then becomes 

2n 



E I 
s 

2R^ 



2it .2 

n o u 



P E A 

/ s s 

* 2R 



0 



Sv 2 



u 



(2-54) 
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The kinetic energy of the shell can be expressed as 



2n 




where 

m = mass of shell per unit surface area 



(2-55) 



2.21 Dilatational Wave 

In view of the symmetries of the problem the displacement com- 
ponents of the shell can be expanded in terms of the following Fourier series 



Ug(S,t) 



00 



Z a (t) cos nQ 

sn^ 

n=0 



V 

s 



(0,t) 



b (t) sin nS 
sn' 



n=l 



(2-56) 



Equations for the strain and kinetic energies of the shell are then expressed 
as the following quadratic forms in the generalized displacements and 
velocities 



U = 



- 



E lit E Art 



E Irt 



'• 5 



) a^ + 



R ' so 



2R- 



2 , 2,2 

) + -sr 



00 

/ (nb + a )' 
/ , sn sn' 



ri=l 



n=l 
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^ t - o 2 , mJtR 

T = m:itRa + — 
so 2 



/'»2 »2 X 

(a + b ) 
sn sn 



n=l 



The Lagrangian equations of motion in terms of the generalized displacement 



coefficients a and b are 
sn sn 



d / ^T \ SU 

dt ^ ^a 

sn sn 



aa ■ ^n 



d ,aT 



au - 



dt ^aa ^ ab *^n 

sn sn 



(2-58) 
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where and are generalized forces corresponding to the displacement 

terms a and b , respectively. The following equations of motion result 
sn sn^ 



E I E A Q 

s s X _ o 

^so ^ ^ t, 4 ''' _2' ^so ” 2rtmR 

mR mE 



El 2 E A 

a + (l-n ) a + (a + nb ) 

sn _4 ^ ' sn „2 '■ sn sn' 

mR mR 



Q 

n 

mirR 



nE A Q 

b ■ — — (a + nb ) = 

sn Jd sn sn mitR 

mR 



n = 0 

( 2 - 39 ) 

n = 2^ ... 

n = 2^ ... 



Now to determine the generalized forces and the principle 

of virtual work is applied. Consider a virtual disp.lacement corresponding 

to an increment 6a of coordinate a 

sn sn 

Virtual work, by definition = Q “6a 

^ ^ n sn 

The external forces acting on the shell are 



00 



a 

rr 



Z A cos nQ 
n 

n=0 



a 



vQ 



00 

^ sin nO 
n=l 



where A and B represent the sum of the forces on the shell due to the 
n n 

incoming and outgoing stress waves 



A = A (in) + A (out) 

B = B (in) + B (out) 

The coefficients A^ and B^ have been discussed previously, 
a virtual change of the generalized coordinate a^^ is 

2jt 



Virtual Work = 




“ 2 

A cos n0 • 6a “ 'Rd9 
n sn 



The work due to 
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The generalized force term is determined hy equating the two expressions of 
work, from which we get 

Q = 23tR A n = 0 

o o 

= xR n = 1, 2, . , . 

and hy a similar computation 



Q = JiR B 
n n 



The equations of motion now heccme 



E I 

.. s 

3 , + 1 

sn 

mR 



2 E A A 

(l-n^) a + — ^ (a + nb ) = — 
' sn t^ 2 sn sn m 

mR 



nE A B 

s / _ \ n 

b + — ^ (a + nb } = — 
sn ^2 sn sn' m 
mR 



(2-40) 



2.22 Shear Wave 

Because of the difference in symmetry associated with the shear 
wave as compared to the dilatational wave^ displacement canponents in the 
case of the incoming shear wave are written as 






y a (t) 

/ f sn^ ' 



n=l 



sin nB 



n=0 



(2-ia) 



Expressions for the strain and kinetic energies in terms of the generalized 



coordinates a and b now become 
sn sn 



U = 



E In 
s 

2R' 



00 — 00 

V 2 ,, 2>2 V” V t 



nb + a )‘ 
sn sn 



n=l 



n=l 



(2-42) 
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, = . H|£ X (iL - \n) 

n=l 



The equations of motion determined in a manner similar to that of the 
dilatational wave are then 



a + — 7- (l-n ) 
“ me’* 



E A A 

a + — ^ (a - nlD } = — 
sn ^2 sn sn' m 
mn 



b 

sn 



nE A 
s 



mE 



B 

(a - nb ) = — 
^ sn sn' m 



( 2 - 1 + 5 ) 



2e25 Effect of Additional Mass 

The equations of motion deriveci thus far assumed no mass within the 
shell. The simplest way to consider the effect of additional mass is to 
assume that it is distributed within the shell symmetrically about the axis. 
The total quantity can be assumed equal to 23iRm* ^ with m* the magnitude of 
the added mass measured in terms of the surface area of the shell. If it is 
assumed to move with the same velocity as the mass center of the shell, its 
kinetic energy can be expressed as 

2n 2 



T* = 



2JtRm* 






/ (-u^ cos 0 + V sin G)dQ 



( 2 - 44 ) 



for the case of the incoming dilatational wave. If u and v are replaced by 

s s 

their Fourier expansions, and if the indicated integration is carried out, 
the kinetic energy T' due to the additional mass may he represented as 



T' = 



m' itR 



^%1 



- 2a b 
sl si 






This kinetic energy must be added to the kinetic energy of the shell which 
was derived earlier. Since the above expression affects only the n = 1 mode, 
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the other modes need not he considered. With this additional kinetic energy, 
the equations of motion for the n = 1 mode become 



E A 



m 



^sl 2m ^^sl ‘ ^sl^ „2 ^^sl ^sl^ m 

rrm 



m' 



E A 
s 






^sl ■ 2m ^^sl ■ ^sl^ ^ „2 ^^sl ^ ^sl^ m 

miA 
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If now we add the above equations, we get 



2E A 






(a + b ) = -;;r + 



^sl ^sl ^2 ^'"sl ' “si' m ■ m 
mK 



or 



a’ = “^ n 
si si 



2E A A, B, 

s / \ 1.1 

^2 ^ si si' m m 

mK 



and 



b , = - a 
si si 



2E A A, B, 

s / , \ 1 1 

^ (a + b ) + ““ + 

^2 si si m m 
mR 



By making the appropriate substitutions for a^^ and 'the equations of 

motion can be written as 



E A 



B 



®sl ^2 '“si ' “si 
mR 



1 



m 



s / V, i 2m+m' 

“ m ■ 2(m+m' ) m 2(m+m' ) 



(2-46) 



E A 
s 



B, 



^ s / -u ) _ ^ ^ 2m+m* 

°sl ^2 ''^sl “si' “ m ‘ ^ra+nt7 m ' 2(m+m' ) 



Note that these equations are similar in form to the equations derived earlier 
without the additional mass; the effect of the additional mass merely alters 
the right-hand side. 

In the case of the incoming shear wave^ the kinetic energy of the 



additional mass is 
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2n 



T' = 



2nRm’ 1 
2 



/ (li sin 0 + V cos d) dS 

-^0 " " J 



(2-47) 



and by a similar method, the equations of motion for the n 
written 

E A 



= 1 mode are 



A B 

_i . ^+m' _ _1 . m' 

si ' '“si “si' ^ m 2(m+m') m 2(m+m' ) 

mR 



a , + — ^ (a - b ) - 



A. 



m' 



B, 



(2-48) 



2m+m' 



EsA 

^sl ^ ^^sl " ^sl^ "" " lir ’ 2(m+m' ) "in’ ' 2(m+m'' ) 
mR 



The significance of the additional mass on the numerical results is discussed 
in a later chapter. Inclusion of a flexible support for the additional mass 
is also possible by changing the above equations; however, this problem is 
not treated here. 



2. 3 Boundary Conditions 

It is assumed that the shell is attached to the medium permitting 
no differential displacements at the boundary between the two. Thus, con- 
tinuity of stresses and displacements are maintained at the boundary. A 
convenient way of satisfying the boundary condition of equal displacements 
of the shell and medium is to equate the corresponding velocity components 
of the shell and medium. This is done modewise for both the dilatational 
and shear incoming waves by merely equating the coefficients of the velocity 
terms as follows 

a = a (in) + a (out) 

b = b (in) + b (out) 

sn n ' n ' 



(2-49) 
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2 . 4 Summary of Equations in Non-Dimensionalized Form 

It is convenient to express the equations derived earlier in 
non-dimensionalized fonn by introducing certain dimensionless variables and 
parameters . Let 



a 

n 



a 

sn . 
R 




P 

n 



b 

sn 

R 




T 



tc. 



R 



(2-50) 



be the non-dimensionalized form of the displacement and time variables. 
Parameters are expressed as 



E 

s 

E 



P 



(l+v)(l-2v) 

— — 



R 

A 

R 

I 



R 



5 



( 2 - 51 ) 



The last two parameters given above become unnecessary when considering an 
unstiffened shell since for this case A and I are given by 



A = t and 



1 ^2 ^ 



(2-52) 



Although the equations were derived considering an unstiffened shell of 
uniform thickness^ the response of a shell with thin, closely spaced stiffeners 
can also be approximated by specifying A and I separately from t. The 
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equations of interest in non-dimens ionalized form are 

d = PI (0, ) + Q1 (F,G) 
n n' 1' n' ■’ 

^n = Q2^(F,G) 

(2-55) 

a = Nl (a ) + P5 (9,) + Q3 (F,G) 
n n' n^'^n'^ n' 1' -^n' ■’ 

p = N 2 (a ,p ) + Pi^ (0j + (F,G) 



where PI , P2 , P5 , and Pi+ are functions whose values are determined 

directly from the position of the incident wave. *^^n^ ^^n^ *^^n 

functions of the outgoing shear and dilatational waves: and N1 and N2 are 

n n 

functions of the displacements. These various functions in non-dimensionalized 
form are derived from equations given earlier. Given here are those for the 
case of the incoming dilatational stress wave. 
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2(l+v) sin n0 



n = 0 



n = 2 
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2 • 5 Other Equations of Interest 



Dilatational Wave 

The hoop stress in the shell is detemined from the relation 



a 



GQ 




e 



Q 



where 



e 



- 1 (__£ + 
a ■ R ^ 



u ) 

s 



By substitution of the Fourier expansion of the displacement components the 
hoop stress can be expressed modewise as a fraction of the absolute amplitude 
of the incoming stress wave, as 




I 




+ a 
n 



cos nG 
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The bending stress is determined from the relation 

E d 

S / S . \ 

bG^ 

which may also be expressed modewise as 





(n^-l) a 



cos nG 



n=0 



( 2 - 56 ) 



where d is the distance from the neutral axis of the shell to its extreme 
fiber o 



Stresses within the medium at any radius r may be determined by 

summation of modal stresses as follows 

00 ^ 
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( 2 - 57 ) 
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00 
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i"pi 
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are 
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sin nQ 



n 



cos nQ 



n' n n " 

stresses of the incoming and outgoing stress waves. 



Shear Wave 

The corresponding equations for the incoming shear wave are^ for 

the shell 




(a -np ) sin n9 
V n 



n=l 



d 2^ (l-n ) sin nO 
n=l 
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and for the stresses in the medium 
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sin nQ 
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00 — 

B 



n 
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sin nQ 
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where all terms are as defined "before. 
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CMPTER III 
METHOD OF SOLUTION 



5.1 General 

In this chapter are presented the numerical techniques used in 
solving the problem. Since the methods employed are similar in theory for 
both the incoming dilatational and shear waves ^ only the solution to the 
incoming dilatational (compressional) wave will be discussed in detail here. 
The equations (2-55) "to be solved consist of two pairs of coupled integro- 
differential equations in the generalized coordinates of the shell and the 
potential functions. The integral teirnis in these equations contain elements 
of the outgoing wave potentials^ both dilatational and shear ^ whose values 
must be determined through application of the boundary equations at each 
instant of time considered. 

In the numerical solution^ time is taken to be the independent 
variable^ and is expressed in terms of the half transit time of the incoming 
wave across the cavity, t = tc^/R. Increments of time are expressed as a 
fraction of the half transit time; if N is defined as the mmiber of time 
steps for the wave to travel one radius, then i/N denotes the elaspsed time 
after the ith step has been taken (Fig. ^4-) . 

The position angle, 0^^, can be defined for any step in time as 

0^ = arc cos (l - i/N) (3“l) 

After total envelopment, i.e., for i/N > 2, 0^ becomes equal to jf. 

5‘2 Numerical Integration of the Potential Functions 

Lamb (^4-) notes that if a point source f(t)dz is located at the 
point (x=0, y =0, z =z), its effect at a distance r from the origin in 
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the xy plane can be represented by the equation 

n I 2 ^ 2 

cp = — a — f(t - ^aut2_) az 

% r^+z^ 

which may be represented graphically as in Fig. 5- To obtain the effect of 
a line source of density F(t) on the z axis^ spherically symmetric point 
sources with the same variation in time can be assumed to be situated all 
along the axis^ and their combined effect represented as an integral over all 
the point sources 



1 








(3-2) 



The limits of minus infinity to plus infinity are shown here in order to 
represent the general case. It will be shown that for a disturbance beginning 
at some definite time^ the limits of the integral may be taken as finite. 

Now let 

2 2 2,. 2, 2x 

r + z =r(.l + z/r; 

2 v,2 

= r cosh u^ 

and the above integral can be written as 

00 

9 = / F(t - cosh u ) du (3-5) 

J 0 J. j. 

from which the graphical interpretation in Fig. 6 is made. 

Consider for the present the determination of the integral values 
at the boundary r = R. We define t = 0 to be that time when the first effects 
of the outgoing stress waves (i.e.^ from the point source at z = O) reach the 
boundary. As time increases the effects of other point sources progressively 
further away along the z axis will reach the boundary in the xy plane. Thus 
it can be seen that though the limits of infinity are given in the integrals^ 
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for finite times the limits can he represented by finite values. The maximum 

value of the variable R cosh u^ in equation (3-3) that need be considered can 

therefore be represented as a function of time. At t = 0 

R cosh u, = R 
Im 

that is, u, is zero for time zero. For any time later 

^ I m 



R cosh u, = R + tc^ 
Im 1 



in the case of the outgoing dilatational wave (Fig. 6). R cosh locates 
the point farthest along the z axis whose wave front has just reached the 
boundary in the xy plane^ c^ being the velocity of wave propagation. In 
non-dimensionalized form this can be written 

cosh u^ =1+T=1+^ 

Im N 

where the terms are as previously defined. Since the integration is over all 
point sources whose wave front has reached or passed the boundary in the xy 
plane^ the integral expression of equation ( 5 - 3 ) can be written with new 
limits 
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Im 

= / F(t - — cosh u. ) du^ 

^0 "1 ^ ^ 



Let 



cosh u^ = 1 + ^ 



(5-4) 



where ^ now becomes the variable of integration varying in value from 
0 to T. In terms of ^ the integral expression becomes 



cp = 



\ F(-l+T-0d^ 

0 n / 1+^/2 



(3-5) 



Note that a singularity of the integrand occurs at ^ = 0^ for all values of x. 
Paul (6) has presented a numerical technique for evaluating integrals of this 
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general type taking into account the singular point. The method is essentially 
a modified trapezoidal rule in which weighting coefficients are obtained. If 
the range (O^t) of the integral is subdivided into subranges each of equal 
length At^ the above integral for any time t = IAt can be represented as 



cp = 




m=0 




+1)At 



mAx 



\/2^ 



( 3 - 6 ) 



where f(0 is a polynomial of first degree in ^ which approximates the 
quotient 

F(-l+t-0 

n / 1+^/2 

as a linear function over each subrange. The results of the integration of 
equation (5”6) are then put in the form 

(i-i) 
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AM • R • F. + M • R • F. 
m m i-m m m+1 i-m+1 



( 5 - 7 ) 
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J L 
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represent a set of weighting functions which can be computed for any interval 
being considered. For large values of m these functions were written in a 
form more convenient for computational purposes by asymptotic expansion^ to 
get the following result 
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For an illustration of the variation of AM and BM functions with see 
Figo 10, is termed the multiplying factor which for this particular 
integral is equal to 




All of the various integral expressions for the potential functions were 
simplified to the form of equation (5‘’6) thereby permitting use of the above 
weighting factors for all integrals. However, the expression for the multi- 
plying factor R^ varies for the different integrals. 

The integration for the outgoing shear potential t is handled in a 
similar manner, except for a minor modification which results from the fact 
that the shear wave travels outward with velocity c^. 

Consider for illustrative purposes some shear potential represented 
as 

00 

^ = / c(t - — cosh u^) cosh u^ du^ 

Jq Cg 2 2 2 

Since the incident wave strikes the shell at time t = 0, the shear wave front 
must also reach the boundary at this instant in order to satisfy boundary 
conditions. Thus at t = 0 

R cosh u^ = R 
2m 

and for any time later 

R cosh u_ = R + tco 
2m ^ 

using the same reasoning as for the outgoing dilatational wave. We can express 
the above equation in non-dimensional ized form as . 

cosh u^ = 1 + k T 
2m c 
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where 

O 

is a constant relating the velocities. 

Again since is a variable of integration we can write 

cosh = 1 + k ^ 

2 c 

where ^ as in the case of the dilatational wave is the new variable of 
integration with limits 0 to T . In terms of this variable the integral is 
now written as 






r G(- i + T - 0 



0 






1 + 






(5-10) 



Note that the shear wave starts from the z axis at time T =-l/k which is 
earlier than for the dilatational wave since it travels at a slower speed 
(with fictitious material within the shell) and must reach the boundary at 
time equal to zero. The numerical integration of equation (j-lO) is 
accomplished in a similar method to that which was done for the dilatational 
potential. The result can be expressed as 

(i-i) 



- I 

m=0 




G. + BM 
i-m m 



■ Sti+1 ■ 



0. 

1 -m+1 



(5-11) 



where AM^ and are weighting factors identical to those derived previously^ 
and 0^ for this particular integral is 






n/V ' (l + k ^ ) 



cm 



1 + 



k ^ 
c m 



The expression for ft is again dependent on the form of each integral considered. 
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For an illustration of the variation of the F and G functions with time^ 
see Figs. 11 and 12. 



5-5 Solution of the Basic Equations 

The equations of motion for the shell are solved by an iterative 
procedure known as the Newmark Beta Method ( 5 ) with which values for the 
potential functions^ and the accelerations^ velocities^ and displacements 
of the shell are determined. 

In general the method consists of using a step-by-step integration 
technique over successive time intervals assuming a specific variation of 
acceleration during each interval. If we assume that at time t = (i-l)AT all 
values of potential functions^ accelerations^ velocities and displacements 
are known, the method becomes that of determining the corresponding values 
at a time t = iAr . For an assumed linear variation of acceleration over each 
interval the equations of interest are 



a . = a 

n,i 



n, i-1 






n, i-1 



+ q: . ) 



n,i' 



• At 

a . - cc . ^ + Ata . ^ + -7- (2a . ^ + a . ) 
n,i n,i-l n,i-l 6 n,i-l 



* • At /** •" \ 



(3-12) 



At 



P . = P . + ATP . + “7 . n 

^n,i ‘^n,i-l ‘^n,i-l 6 ^ ^n,i-l 



(2p . T + p .) 



where n denotes the mode and i the time. Note that in the above equations 

the velocities and displacements are expressed in terms of a . _ and 6 . _ , 

n,i-l *^n,i-l^ 

acceleration components which generally are unknown. 

Values for the acceleration components at time t = iAT are assumed 
and the velocities and displacements determined from the above equations; then 
by use of the boundary equations to determine the values of the potential 
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functions^ values of the acceleration components can be computed from the 
equations of motion. The computed values are then compared with the assumed 
values and if they agree or are within an arbitrarily established tolerance^ 
the step is completed; if not^ the values of ^ and ^ just computed are 
taken as new assumed values and the cycle repeated until the criterion is 
satisfied. A detailed discussion of the procedure follows. 

The equations of continuity at the boundary at time t = iAr are 
given in equations (2-49) rewritten here in the form 



a . = PI (aj + 
n^i n^ 1 ' 



iAt 

f 
0 

IAt 






( 3 - 15 ) 



P . = P2 (0-,) + / Q2 . .F.; G„,G-, , . . .G. ) 

^n^i n^ 1' J n^ 0^ 1^ 0^ 1^ ^ 



0 



where PI and P2 denote functions related to the particle velocities in the 
n n 

medium at the boundary due to the incoming stress wave; their values may be 
determined directly at any time. The integral expressions denote the con- 
tribution due to the outgoing waves . All values of F and G are known except 
F^ and which are to be determined from the above equations. These equations 
are now written isolating the unknown values of F and G terms 

IAt 



6c . 

n,i 



= PI (0, ) 

n 1 









P. . = 



At 



At 

+ / Q1 (F. ^,F.; G. ,,G.) d^ 

Jq n^ 1-1^ 1 1-1^ i' 

IAt 

At 
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/ Q2 (F. ,,F.; G. ,,G.) d^ 
Jq n' 1-1^ 1-1^ 1' ^ 



The integrals with limits At to iAT are evaluated using the numerical 

technique described in a previous section. 

The next step is to assume values for a . and 6 . , the radial and 

^ n^i n^i^ 

tangential components of acceleration of the shelly from which velocity and 
displacement components are determined using equations (5»12). With these 
values known^ equations (j-l^) become effectively two equations with the two 
unknowns^ and which can then be determined. 

The equations of motion are given in equations (2~55) here 

written in the form 



iAx 



a . = Nl(a p .) + Pj (0j + 
n,i ^ n,i 1 ' 



6 . = N2(ci p .) + P4 (0j + 

n,i n' 1' 



^^n 0^1^ 1*^ 0“^!^ i' ^ 



(3-15) 



iAT 



Q4 (F^,F^ , . . .F.; G^,G, ,...G.) 

n o^ 1^ 1^ o'^ 1^ i' ^ 



from which the acceleration components oc . and p . are now computed. If 



n,i 



n,i 



the computed values agree or are within an arbitrarily established limit of 

the assumed values^ the step is completed; otherwise the computed values are 

used as the assumed values for the next cycle of iteration. 

The procedure described above requires that all parameters at t 0 

be known, including the initial values F and G » Initial values of velocity 

o o 

and displacement of the shell are specified to be zero. From the short-time 
approximation presented in Section 5*7 initial values of acceleration and G^ 
were determined to also be zero, and the values of F^ to be the following 



F = 
o 
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TT 



-1 



2(-l) 



n+1 



n = 0 
n = 1,2, 



(3-16) 



The short-time approximation also indicated that the radial acceleration 

1/2 ^/2 

components near t = 0 varied as t ' and the tangential components as t ' ; 



^5 



thus to improve the accuracy of the calculations^ equations (5“12) for the 
first step in time were modified to 



^ 1 
n,l 
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n,l 
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Without this modification^ the results for the first step cannot be brought 
into acceptable agreement with the short-time approximation o 

The n = 0 mode is simplified somewhat because it has no tangential 
component of displacement^ and contains only the dilatational component of 
the potential functions . 



5 • ^ Stresses of the Shell and Medium 

Solution of the equations of motion yields modal values of the 
functions F and Gy and the modal acceleration^ velocity and displacement 
components of the shell. The stresses in the shell are determined directly 
from the displacements using equations (2-55) £^^d ( 2 - 56 ). 

Stresses in the medium at any radius and time are determined using 
equations (2-5?) • They are here written in the form 



g(r^ iAr) 
a 

P 



Z -X 



T (F)d^ 

nx / ^ 



W^(G) 



n=0 



( 3 - 18 ) 



where represents the stress due to the incoming wave and the integral 

terms that due to the outgoing waves . All are functions of both radius 
and time. Consider the stresses in the medium for any radius equal to r at 
time T = iAT = i/N. 

9^ (Fig. 7) is now determined from the equation 
0^ = arc cos — (l - i/N) 



(5-19) 
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and takes on values from 



^ =- 
1 r 



for T = 0 



to 



= TT 



for T > (1 + i) 



The integral terns must be recomputed for each radius and time considered. 
As an example consider the integral 

00 

cp = / F(t - — cosh U-, ) du-, 

J 0 -‘- 

and let t = ~ - 1, where k = R/r, represent the nondimens ionalized time 
K. r 

r 

after t = 0 required for the outgoing dilatational wave front to reach the 

radius r being considered. Then for t < the integral is equal to zero^ 

and for t > t the maximum value of fr cosh u., ) can be written 
r 1' 

r cosh U-, = r + t ’ C-, 

Im 1 



or 



cosh u^ = 1 + k T ' 

Im r 

where t' = represents the time after the wave front has reached the 

radius r. We express cosh u^ as 

cosh u_ = 1 + k C 
1 r^ 

where ^ is the same variable of integration as was used previously. The 
integral is then expressed in terms of ^ as 



cp = 




F(-1+t'-0 



( 3 - 20 ) 



The numerical method of integration used is identical to that described 



^5 



earlier for the integral at the boundary^ with the result that the same 
weighting factors AM^ and are applicable. Note, however, that the limits 
of the integral are now 0 to t*; also the multiplying factor for this 
particular integral becomes 



R 

m 





r m 
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By a similar analysis, the integral of the shear potential 
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= / G(t - — cosh u ) cosh Up du 

J A a d. 






can he put in the form 
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where 



T ' = T - T 
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For the potential t 
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The weighting factors AM^ and BM^ remain the same but the multiplying factor 



0^ for this particular integral becomes 
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5-5 Time Dependent Stress Wave 

The problem considered thus far has dealt with an incoming stress 
wave with a step distribution in time. The results obtained for the step 
wave can be used through the application of Duhamel*s integral to find values 
for the response of a shell and medium to incoming waves with any time 
variation. As an example^ a stress wave which decays exponentially with time 
(Fig. 8) according to the following equation is considered 

-k ^ 

a (t) = a (l - e ^ (5-22) 

p\ / pq\ / 



represents the time at which the stress wave decays to zero and k is a 
parameter which is related to the shape of the curve. Stresses at any time 
equal to iAt can be found by the application of Duhamel*s integral^ here 
written in the form 



a ( iAt ) = a ( iAx ) + 

Su 




da (t) 
P 
dT 



a (iAT-t) dT 
sx> 



(5-25) 



a is defined to be that stress resulting from an incoming step wave of 
Su 

amplitude in equation (5-22). In the numerical analysis^ the time 
dependent wave is approximated by a series of rectangular sections as illu- 
strated in Fig. 8. 



5*6 Description of the Computer Program 

The problem was programmed modewise for a high speed digital 
computer (CDC l604) using Fortran language^ and considering only the first 
three modes. Input data consist of the following parameters: 

(l) The time intervals at which computations are to be performed^ 
expressed as the number of intervals required for the incoming wave to travel 
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one radius (one-half transit time). The degree of accuracy achieved is 
dependent on this parameter^ the smaller the interval the more accurate 
the results. However^ the machine time required for a given mmiber of transits 
of the incoming wave varies inversely as the square of the interval size ap- 
proximately^ thus some sacrifice in accuracy is necessary to reduce the time 
of computations required. 

(2) The total time over which the computations are to be performed^ 

expressed as the total number of time intervals to be considered. Generally 

speaking^ all values seemed to have reached their asymptotic (static) values 

within ten transit times of the incoming wave across the cavity. 

(5) The ratio of the moduli of elasticity, E /E, where E is the 

s s 

“plane strain** modulus for the shell and E the modulus of elasticity of the 
mediimi . 

(k) The mass ratio, p /p, where p is the mass of the shell and p 

s s 

the mass of the medium per unit volume. 

(5) Poisson’s ratio of the medium. 

(6) The ratio of the thickness of the shell to its radius. When 
considering a shell whose a.rea. A, and moment of inertia, I, are not directly 
related to the thickness, A and I must be specified separately. 

(7) The amount of additional mass within the shell expressed as a 
fraction of the mass per unit surface area of the shell. 

( 8 ) The number of radii to which stresses in the medium are 

desired. 

(9) The time intervals at which output data are desired. 

(10) The angular increment at which output data is to be computed. 
Because of symmetry only values between 0 and I80 degrees need be considered. 
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Additional input quantities in the case of the exponentially- 
decaying stress wave are and parameters which define the duration and 
shape of the stress pulse ^ respectively. 

Output data consist of the acceleration^ velocity displacement^ 
and stress components of the shell for specified angles and times; and the 
stresses in the medium for specified radii; angles and times. 



5.7 Short-Time Approximation 

As a check on the accuracy of the machine solution for short times ; 
the houndaiy equations and the equations of motion were solved approximately 
hy making a series expansion of all pertinent functions in terms of time as 
the independent variable. Although the following discussion is limited to the 
case of the incoming dilatational wave; the basic principles are the same for 
either of the types of wave considered. 

The expressions for the velocities and stresses in the medi-um 
around the boundary due to the incident wave all have been written thus far 
in terms of 0 ^; the position angle of the wave. In terms of nondimens ionalized 
time T; all functions of can be written in terms of t using the following 
relations 



cos 0 ^ 
sin 



= 1 - T 



= 




= '/ir 



1 + t + 



_ 1 _ 
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( 3 “ 24 ) 



The integral values representing the effects of the outgoing waves 
are also expressed as functions of t. The following example will illustrate 
the technique used to accomplish this. For example; consider the transforma- 



tion of the integral 
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sv 

CD = / F(t - — cosh u. ) du^ 

Jo ^ ^ 



As was shown in Section 5*2; cosh in the integral varies in value for any 
given time r from 1 to 1 + x . It is convenient for purposes of analyzing the 
integral to represent the variable cosh u^ as follows 

cosh u^ = 1 + cjot ( 5-25) 

where t is now a fixed value in the integration and co is defined as the 
variable of integration whose value ranges from 0 to 1. The function 



F(t - — cosh u^ ) = F(-1 + t(1-co)) 

Cl 1 

is expanded in terms of a power series as 

00 

F(-l + t(1-u))) = F + ) 

O / I 1 

i=l 

where 7. are unknown coefficients of the series. Since 
1 



du. = 






1 - -p + (wt) + . . . 



'*52 
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the integral can now be written as 

1 
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= il f 2^ 
J 2 *^0 ''/co 
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^ Cx)T 
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(5-26) 



from which after performing the indicated integration 



cp 



= '/ir 



F + t(| 7, - z F ) + ( 

o 5 1 0 o ' 



(3-27) 



The integrals which represent the effects of the outgoing shear 
wave can also be transformed in somewhat similar manner. Consider for example 



the integral 
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ilr = / G(t - ^ cosh u^) cosh du^ 

Jq ^2 ^ ^ ^ 



In this case cosh Ug varies from 1 to 1 + k^T where 
Therefore we write 



=2 , 

— and T = — 
c^ R 



cosh = 1 + mk T 
2 c 



(3-28) 



where w is now the variable of integration ranging in value from 0 to 1 in 
the integral. The function 



G(t - ^ cosh Ug) = G(- + -rCl-w)) 

^2 c 



is expanded in terms of a power series as 



G(- + tCI-u)) = g^ + ) e^T^(l-u)- 

i=l 



where are unknown coefficients of the series. From equation (j-28) 
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The integral written in terms of the variable u) is now 
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and performing the integration 



1 + uk T 
c 



dm 



♦ ^ “o> * 



(3-29) 



The same basic technique is applied to all the integrals so that now the 
effects of both the incoming and outgoing waves can be written in terms of t , 
Substitution into the continuity equations and the equations of motion yield 
the following for any mode 
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» . 1 

V ^ ? 

a = 2^ fi.(F,G)T 
i=l 

i — 

P = ^ f2^{F,G)i ^ 

i=l 

(3“50) 

00 i 1 

a + il.(a,p) = ^ f3.(F,G)T ^ 

i=l 

p + £2^(a,p) = ^ f4.(F,G)T ^ 

i=l 

lere fl., f2., 3.nd fk. are coefficients which contain certain elements 

1 1 1 1 

7^ and of the potential functions; il and 12 are known functions 
' the displacements. 

To solve the above equations, the displacement components of the 
Lell are expressed as Frobenius (8) type series 



<x> 




i=l 



00 




i=l 



(5-31) 



lere p^, q^, and c are unknown coefficients. Substitution into equations 
)- 30 ) results in the following set of equations 
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(3-32) 
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The coefficient c is now determined by inspection. Then through a step by 
step process which involves the equating of coefficients of like powers of 
values of p^^ q^, 7 ^^ and are determined. The following equations are then 
used to find values of the potential functions^ and the displacement components 
of the shell. 
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G + 7 e.T 
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c+i 



i=l 



(3-53) 



Velocity and acceleration components may be determined by differentiation of 
the above expressions for the displacements. 
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CHAPTER IV 

DISCUSSION OF RESULTS 



General 

Results of computations performed to determine the effect of the 
various parameters are discussed in this chapter. 

Although equations presented throughout the study have been written 
to include an infinite number of modes, the greater part of the actual calcu- 
lations performed and presented here are the results obtained considering only 
the modes n = 0, 1, and 2. It is important to note that during envelopment of 
the shell by the plane stress wave, a Fourier series representation of the 
incoming wave is objectionable in that the series at this stage is slowly con- 
vergent, thus necessitating a large number of modes to accurately represent the 
plane wave. However, after passage of the wave across the cavity, the Fourier 
expansion of the incoming stresses around the boundary results in coefficients 
of all modes except n = 0 and 2 becoming identically equal to zero for the 
plane dilatational wave, and coefficients of all modes except n = 2 becoming 
identically equal to zero for the plane shear wave. Therefore, stresses due 
to the outgoing waves in modes corresponding to those of the incoming wave 
whose coefficients become zero must also eventually vanish at long times. The 
limited study conducted for modes greater than n = 2 indicated that the 
maximum effect of the higher modes occurs within one transit time of the 
incident wave and rapidly decays, thus contributing relatively little to the 
maximum response of the shell which occurs after several transit times. 

However, for determining the early time response of the shell and medium, 
the higher modes are significant and should be considered in further extension 



of this work. 



In the tables and figures to be discussed, quantities given in 

non-dimensionalized units are defined by equations (2-36) and (2-5O) . Stresses 

are given in units of the absolute value of the amplitude (|a | or ja 1 ) of the 

P s 

incident wave; a negative stress means a compressive response to an incoming 
compressional (Fig. l) or a positive shear wave (Fig. 2 ). The physical 
properties of the shell relative to those of the medium, as well as the nature 
of the incoming wave are indicated on the graphs. Unless otherwise stated, 
the shell is considered to be an unstiffened one so that its cross sectional 
area and moment of inertia are related to the thickness as given by equations 
(2-52). Except where indicated, there is assumed to be no additional mass 
within the shell. Numeral subscripts denote the mode number. 

The shell and medium have been assumed to exhibit linearly elastic 
behavior throughout their stress histories, which for the practical problem 
does not permit evaluation of any spalling or non-elastic effects. 

Values of stresses given are in addition to those which exist prior 
to the arrival of the incident wave. For the elastic case, the effect of 
prior stresses such as those resulting from the overburden may be taken into 
account by merely adding them to stresses caused by the incident wave. 

For clarity of presentation and because of the impracticability of 
including solutions for all possible permutations of the parameters involved, 
the discussion in this chapter is limited to a few representative cases. 

kc 2 Modal Response of the Shell and Medium 

Figures 11 and 12 illustrate the shape of the modal components of 
the dilatational and shear potentials obtained in the solution to a typical 
problem. It appears that a singularity occurs at one transit time in the 
case of the F functions resulting in the slight irregularity of the curves at 
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this point. Since computed values of stresses were determined to be rather 
insensitive to relatively large variations in values of the F and G functions ; 
the effect of the irregularity would seem to be slight. 

Figures I 5 through 20 show modal acceleration^ velocity^ displace- 
ment and stress components for the shell and stress components for the medium 
at the boundary^ as they vary with time. Static values shown were computed 
using the method given in Appendix B. 

The high accelerations computed near the beginning are not truly 
representative of the actual case^ since they are the result of assumptions 
made earlier in deriving the equations for the shell. The shell was repre- 
sented by a line describing its middle surface which permits no variation in 
accelerations^ velocities^ and displacements of particles through the actual 
thickness. Also^ no provision was made for refraction of the incident wave 
through the shell lining. These limitations restrict the applicability of 
the solutions to a shell whose thickness is small relative to its radius. 

The n = 1 mode is primarily a translational one which accounts for 
the rigid body translation of the shell after it has been enveloped by the 
incident wave. Thus^ it can be seen that the velocity components for this 
mode approach constant values equal to the velocity of the medium behind the 
incident wave fronts and displacements grow without bound reaching a straight 
line variation with time. Note that the stresses contributed by this mode 
reach their peak values within one transit time and quickly damp out^ 
approaching zero asymptotically. For the incident shear wave^ the n = 0 mode 
is also a rigid body movement which accounts for rigid body rotation^ and 
contributes little to the stresses. 

Modal quantities obtained are coefficients of Fourier series; 



therefore^ the total response or effect is determined by adding the coefficients 
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multiplied by the appropriate sine and cosine terms for any desired angle. 
Figures 21 through 24 show the time variation of stresses in the shell and 
medium for various angles when the first three modes are summed c The maxi- 
m\om stress in the shell for any angle may be determined by adding the bending 
stress to the hoop stress. This is indicated in Figs. 21 and 24 by the 
dotted line above the hoop stress. 

Figure 25 is given to illustrate the relative magnitudes of the 
hoop stresses in the shell and medium for several thicknesses of shell. This 
also shows the effect of varying the relative thickness of the shell on the 
hoop stress in the medium^ The dotted line indicating the hoop stress in 
the mediiom for an unlined cavity was obtained from the report by Paul (6) . 

Figures 26 and 27 show how the relative thickness of the shell 
affects the radial and shear stresses in the medium. 

As was discussed earlier^ stresses in the medium for any radius 
can be determined by reevaluating the integral terms which represent the 
effects of the outgoing waves ^ and adding them to the Fourier expansion of 
the incident wave. Figure 28 shows the modal and total radial^ hoop^ and 
shear stresses which were computed for a time equal to 10 transit times » 

These values are compared later with the static stresses^ but on this figure 
the static stresses do not differ by more than the thickness of the lines ^ 
and therefore are not shown. The time variation of the radial and hoop 
stresses in the medium for various radii are shown in Fig. 29 for the incident 
dilatational wave^ and in Fig. JO for the incident shear wave. 

4. 3 Short-Time and Asymptotic Comparisons 

A method for obtaining a solution to the problem which is accurate 
for veiy short times (t « l) was presented in Section 5«7* This was desirable 
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to validate the machine solution and to determine the effect of the size of 
time interval selected. The results for a representative problem are shown 
in Figs. 51 through 53 the incident dilatational wave^ and in Figs. 
through 56 for the incident shear wave. Four terns of the series representing 
the F and G functions^ and three terms for other quantities were used in the 
short-time solution. 

As can be seen from the graphs^ good agreement was obtained for 
very small values of time^ somewhat shorter time being obtained for the shear 
wave as compared to the dilatational wave. The shorter time results from the 
nature of the forcing functions (Eqs. 2-19 and 2-52) which indicate a more 
rapid rise in the incoming stresses for the incident shear wave. 

Within the range of time for which the short-time solution is valid^ 
decreasing the size of time interval for each step of the machine solution 
resulted in closer agreement between the two methods, as is to be expected. 

It also indicated that the stresses and displacements are not as sensitive to 
variations in the interval size as are the F and G functions. 

At the other end of the time scale, i.e., at a relatively long time 
after passage of the incident wave front across the cavity, another check on 
the accuracy of the machine solution is afforded by the asymptotic approach 
of all values to the static results . Figures presented thus far have shown 
that the static condition is approached well within ten transit times. 

The time interval used in the machine solution affects the stability 
of the results for long times. This is indicated in Fig. 57 which shows the 
variations in computed values of the displacement components for mode n = 2 
at relatively long times, for different time intervals. N represents the 
number of time steps required for the incident wave to travel a distance equal 
to the radius of the opening, and 1/N defines the interval size. Note the 








wim 
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smaller graph which shows the percentage difference between the computed 
values at ten transit times and the static values. 

As the interval is decreased the machine solution at long time 
approaches the static solution more closely. Below a certain size of time 
interval^ there is little difference in the results ^ which indicates asymp- 
totic convergence to the correct solution. For this particular problem^ 

N - 50 seems to be "critical” in that for N < wide variations in computed 
values occur. As the mode number increases^ the "critical” value of N in- 
creases rapidly^ and the requirements of computer storage and calculation 
time become decisive factors which make impractical the stuc3y for long times 
of modes much larger than 2. The interval size selected for all problems 
solved (exclusive of the study to determine the effect of the interval size) 
was set equal to l/kO (N = 40) of the half transit time of the incident wave. 

For the static case^ only the modes n = 0 and 2 yield values other 
than zero. Table 1 compares values of shell stresses and displacement com- 
ponents in these modes obtained from the computer solution to a particular 
problem at a time equal to 10 transit times^ with the static solution. Most 
pairs of values differ by less than one percent. Comparable agreement of 
stresses in the medium at various radii are shown in Table 2. Values in 
Tables 1 and 2 were obtained from the solution to a problem whose parameters 
were: = k.O, = J.O, = .05, and v = .25. 

Although the figures and tables presented above were for a particular 
problem^ the discussion given is applicable to all problems which were solved. 
Changing the physical characteristics of the shell and medium within the range 
of values studied had hardly any effect on the degree to which the long time 
machine solution and the static solution agreed. Also^, the absolute value of 
all quantities which should asymptotically approach zero became less than 
.00005 in each case well before ten transit times. 
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1 ^., ij. Effect of Parameters 

Studies were conducted to determine the effect of each of the 
following parameters : 



T] 



E 

t 



E 

E 

i. 

R 




V = V (medium) 



for an unstiffened shell without additional masso Since it was impractical 
to take into account all permutations of the above parameters^ a basic shell 
where 

\ Tip = 5.0 

T)^ = .05 V = .25 



was considered from which each parameter was separately varied to determine 
its effect on the resulting stresses and displacements e Calculations 
performed were only for the case of the incident dilatational (P) wave. 

Of particular interest was the determination of the maximum dynamic 
stresses and displacements (not including rigid body translation) due to the 
incident stress wave. Tables 5 through 6 compare the maximum values obtained 
in the machine solution with the static solution. **DLF/* termed the dynamic 
load factor^ is defined as the factor by which the displacement or stress 
produced by applied as a static load should be multiplied in order to obtain 
the maximum dynamic value. Figures 58 through 4 l are given to graphically 
illustrate the variations in stresses and displacements^ both static and 
dynamic; in the range of parameters considered. Stresses are given in units 
of |a^| and displacements in non-dimens ionalized units defined by equations 

(2-56) and (2-50). 
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Table 5 and. Fig. 58 illustrate the effect of increasing the 
relative thickness of the shell. The range of values selected, for x\ is 
probably much greater than is practical or justified, by the assumptions of 
the analysis, and. was considered, only to determine the trend of the results. 

As Tj is increased, displacements and hoop stresses of both the shell and 

U 

medium decrease; however, the bending stress of the shell and the radial and 
shear stresses in the medium at the boundary increase. Figure 25 shows the 
time variation of the hoop stress in the shell and medium for various thick- 
nesses of shell, including the case of the unlined cavity. 

Increasing r\^ results in a rapid increase in the displacements and 
stresses of the shell and a much lower rate of attenuation of the hoop stress 
in the medium. See Fig. 59 Table 4 . 

Table 5 und Fig. 40 show that the parameter has no effect on the 
static results but does affect the dynamic response. Increasing increased 
the maximum response in both the medium and shell. 

Figure 4l and Table 6 show that within the range of Poisson's ratio 
for the medium considered, as v increases displacements and stresses in the 
shell decrease, with very little additional reduction of the medium hoop 
stress . 

Additional mass within the shell was shown to contribute relatively 
little to the overall response. This is mainly due to the simplifying assump- 
tions which were made in deriving the equations of motion to account for the 
additional mass. Figure 42 shows the effect of the presence of additional 
mass equal to 407 t times the mass of the shell on the dynamic response of a 
particular shell. Additional mass decreases the displacements and increases 
the dynamic stresses in mode n = 1. 
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ko 5 Response to Time Varying Incident Wave 

The results obtained from a solution to an incident wave with a step 
variation in time was shown to he use.ful by application of Duhamel*s integral in 
determining the response to any time varying stress wave. The case of the 
exponentially decaying wave was considered^ the results of which are 
illustrated in Fig« 45 • Substantial reductions in the maximum stresses can 
be expected as the duration becomes smaller. In Fig. 44 is shown the effect 
of a linear rise in the amplitude of the incident wave. Note that for a 
wave with a linear rise followed by a step variation in time^ very little 
decrease in maximum stresses occurred. 



4.6 Comparison with Previous Work 

Baron (l)^ using a different method of analysis, investigated the 
dynamic response of two shells with different physical characteristics sub» 
jected to an incident plane dilatational wave. Figures 45 through 49 compare 
his results for the modal values of stresses and displacements of the thin 
shell, with results obtained by the method of solution given in this study. 
Similar results were obtained considering his so-called stiff shell. Although 
the shape of the response curves can be considered similar in the two reports, 
the magnitude of the dynamic response in Ref. (l) seems consistently higher 
than in the present report, and the long time results asymptotically approach 
values higher than the static solution. The following table compares the 
maximum modal stresses, for Poisson’s ratio of 0 .25> obtained in the two reports, 
and also shows values for the static case. Stresses are given in units of 
I a [. The bending stresses for the stiffened shell cannot be compared directly 



since d, the distance from the neutral axis of the shell to its extreme fiber, 
is not stated in Ref. (l). 
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Quant . 


Mode 


Thin Shell 


Stiff Shell 


Maximum 


Stresses 


Static 


Maximiom 


Stresses 


static 


Ref. 

( 1 ) 


Current 

Work 


Ref. 

( 1 ) 


Cu.rrent 

Work 


Q 

CD 

CD 


0 


-4.79 


-4.56 


- 4.10 


-4.79 


-4.08 


-5.67 


Q 

CD 

CD 


2 


4.69 


4.11 


3.67 


5.42 


5.20 


2.95 


-u 

sb 


0 


0.043 


0.045 


0.039 


0.55 


4 . 08 d 


3 . 67 d 


a. 

sb 


2 


0.140 


0.137 


0.122 


1.70 


15.29d 


ll.89d 



4.7 Conclusions 

Conclusions drawn from the results of the analysis are: 

(1) The method which has been presented is practical for effective 
computation of the dynamic response of a cylindrical shell embedded in an 
elastic medium when subjected to plane dilatational or shear waves. The 
solution presented herein is believed correct since it was checked by independent 
methods at short and long times. 

(2) Peak stresses and displacements in both the medium and shell 
occur sometime after the transit of the incident wave across the cavity; of 
the problems solved^ both the average and mean time at which the peak values 
occurred was equal to 5 transit times. The dynamic effect measured in terms 
of the ratio of the maximum stresses and displacements to the static values 
varied within a relatively small range. The average value of this ratio was 
I0I2 and the mean^ 1 . 11 . Thus^ for the cases considered at leasts the maximum 
stresses and displacements to be expected for a particular situation can be 
roughly approximated by determining the static values and by multiplying them 



by a factor of 1 . 1 . 



(3) The largest stresses in the shell occurred for a relatively 
thin liner in a medium with a low modulus of elasticity^ low mass density^ 
and low Poisson’s ratio. The greatest reduction of the hoop stress in the 
mediim as compared to the unlined cavity results from a relatively thick 
liner in a medium with low modulus of elasticity and high mass density. 
Additional mass within the shell has relatively small effect on the dynamic 
stresses . 

{k) The practical value of tunnel linings to reduce the maximum 
stresses in the medi-um depends on the several conditions mentioned above ^ and 
on the magnitude of the incident stress wave. Under favorable conditions ^ 
reduction of stresses on the order of 50 percent or more is possible. However^ 
for materials such as granite^ smaller reductions can be expected for steel 
liners of practical size. 

( 5 ) Certain assumptions made concerning the behavior of the shell 
have limited the applicability of the analysis to relatively thin liners. 

Future studies of the behavior of thick shells would indicate the effect of 
the approximations used herein. The analysis has also been based on the 
assumption that the behavior of both the liner and medium is linearly elastic 
at all times. Perhaps a more desirable condition would be one in which some 
inelastic behavior is permitted to take place in the medium surrounding the 
shell or in the shell, or one in which some inelastic energy absorbing medium 
such as cinders or foamed plastic surrounds the shell. As a subject of future 
study, it is recommended that the behavior of thick shells, and thin shells 
surrounded by some energy absorbing layer, be considered. 
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APPENDIX A 

DERIVATION OF THE EXPRESSIONS INVOLVING THE POTENTIAL FUNCTIONS 

A .1 General Form of the Potential Function 

The dilatational potential function must satisfy the wave 

equation 



^,2 r^r ^2 ^,2 -^2 ^^2 



(A-1) 



cp in the case of the incoming dilatational wave is an even function of 0 and 
thus can be expressed as a cosine Fourier series 



9 



QO 



(A-2) 



n=0 



where f^ is a function of r and t representing the modal coefficient of the 
potential function. By substitution of the above expression into equation 
(A-l) we obtain the equation that must be satisfied by f . 



n 1 n n 1 n 

^2 r ^r “ 2 ■ 2 :,^2 

or r c^ ot 



(A -5) 



The general solution of f is assumed to be of the form 



f = r^ R 
n n^ ^ ^ 



(A-4) 



Then is some function of r and t satisfying 



2 






■^r 



_il + (2n+l) 

2 r Sr 



(A-5) 



If R^ is a solution it can be shown that the corresponding equation for 



R 



is satisfied by 



R 



n+1 



, Sr 

1 n 

r Sr 



(A-6) 
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By repeated application of this result it can be shown that equation (A- 5 ) 
is satisfied by 



.1 6 



n 



R. = (7 R^ 



n 



where R is the solution of 
o 



(A-7) 




rS^R . 6 r -I 

O 1 O 

^2 r Sr 

-or 



(A- 8 ) 



The solution to this equation for the case of a wave diverging from a center 
is 



R 

o 




r 




cosh u-| ) du^ 



(A -9) 



For the proof of this^ see Lamb Therefore^ the coefficient f^ of the 

potential function is written as 



f 

n 



n/1 S \ 
= r (— 3 -) 

\r 



n 



00 




- — cosh u^ ) du^ 
c^ 1 ' 1 



(A- 10 ) 



For modes 0^ 1;, and 2 ., this can be written 



00 




00 00 

f =— / F"(| ) cosh^u du + ■— / F'(| ) cosh u du 

0 1 '^0 

where primes indicate the derivative of the function F with respect to its 
argument. Note that the expression for the n = 2 mode contains derivatives 
of two different orders. It is convenient to express the function f^ in 
terms of derivatives of a single order which may be done through integration 
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by parts* For the n = 2 mode we then wish to change the second integral to 
an integral involving the second derivative of the function F. Integrating 
by parts gives 



u. 



Im 



^dTl = |•T1 






u 



Im 



T]*d| 



where 



Thus 



I. = F' (e^) 



d| = - ~ F"(|^) sinh du^ 



T] = sinh u 



dT] = cosh u^ du^ 



/ F’ (I ) cosh u du = — / F”(| ) sinh u du + 

o'l 1 11 



u. 



F' sinh u^ 



Im 



The second term can be shown to equal zero since for its lower limit 
sinh u^ = 0; and for its upper limits F' is the integral of the function F^ 
at the wave front which under the assumed initial conditions^ does not exist* 
Therefore, the coefficient of the potential function for the n = 2 mode can 
be written as 



f ^ / F"(t “ — cosh u. ) cosh 2u, du, 

2 ^2 Jo "l ^ 1 



(A-11) 



By a similar process the coefficients of the potential function for any other 
mode can be reduced to the fom 



f 

n n 



n 



F^(t 



cosh u, ) cosh nu, du 



1 1 



(A -12) 



The shear potential function must satisfy the wave equation 
" r ^ \2 %2 J,2 



(A -13) 
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t in the case of the incoming dilatational wave must he an odd function of 9 

and can he represented as 

00 

t = ^ g^(r,t) sin n0 
n=l 

and proceeding exactly as in the case of the dilatational potential the 
general expression for the coefficient g^ is found to he 

00 

g = f G*^(t - — cosh u^) cosh nuo d.u^ (A-l4) 

J Q ^2 ^ 2 2 



A. 2 Velocity Terms 

The equations for the velocity components of a particle on the 
■boundary due to the outgoing waves are 



• - i 

^ ~ Srot ^ r SOot 

1 O Cp d t 
^ ~ r ^0ot ^rSt 



(A -15) 



The components of velocity and the potential functions are expanded in series 
as 



00 

u = ^ a^(r^t) cos n0 



n=0 



V = / h^(r^t) sin n9 

n=l 






00 



n=0 



t = / g^(r^t) sin nO 

n=0 



(A-16) 
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Substitution into equations ^A-15J gives 






• ^ II II II 

^n ~ 6r6t ^ r 5t 

. 6f 6^g 

\(out) = - 7 ^ ^ 



(A-17) 



The expressions for the coefficients f^ and are given in equations (A- 12 ) 
and (A-ll4-)<, Using these in the above equations and by the application of 
integration by parts we get the following velocity terms 



-1^ 



n 



an(out) = - J cosh u^ cosh nu^ du^ 

Cl 0 



111 
' n+1 ' 



+ sinh Ug sinh nu^ du^ 



(A-18) 



b (out) = - i-1- r P^''’^(e ) sinh u, sinh nu-, du^ 
n ' n+1 J ' 1 • 1 11 

Cl 0 



1 ° 



Ac 3 Stress Terms 

Stresses in the medium due to the outgoing waves are written in 
series form as 



rr 



A (r.t) cos nO 



n=0 



TQ 



'QG 



00 

= ^ B^(r,t) sin nG 



(A-19) 



n=l 



00 

= ^ cos n0 



n=0 
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By substitution into equations (2~i0;' the coefficients and are 

expressed in terms of the coefficients of the potential functions, f^ and 



A^(out) = (\ f 2jj.} 



rS^f Sg ng -] 

n n n n 

^ 7 77" ~ '~2 

t-dr r ■ 



rBf 2 Bg 

n n n r 

77 r ®n r n Sr 



B (out) 



o o ^ S 2 1 Sg 1 

2n n 2n °n n 1 n 

rSr^2n -,2 2®n^r Sr 

r Sr r 



(A-20) 



C (out) = 



C^f 






n n ^ ^ n 

f- g - f - n 
or r r n or 



+ \ 



rS f bg ng - 

n n °n 

^ 2 ^ r Sr 2 

-or r * 



Again;, using the equations for anc3 g^ given earlier we can write the above 
coefficients finally as 



A^(out) = -1^^ [' cosh nu. 



n+2 ^ _ 



^ o v,2 

“ + 2 cosh u, 

1. 



du. 



=2 ° 



B^(out) = f sinh 2u^ sinh nu^ du^ 

c^^ 0 



(-l)j 



,n+2. 



n.2 j. ^ '^^2^ >"^2 dUg 

C 2 0 



(A -21) 



Cn(out) = J'^ cosh nu 3 _ 



\ 2 
— ~ 2 sinh Ut 
u 1 



du, 



1 



Li&- /'• 



+ J ^ '^ 2 ^ ^2 ™2 ^^^2 

^2 ° 
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APPENDIX B 
STATIC SOLUTION 



B.l Dilatational Wave 

The static solution presented here is Based on the application of 
the theory of stress functions presented hy Timoshenko ( 9 )- Under static 
conditions^ it can be assumed that at large distances from the boundary of 
the cavity the state of stress in the medium is equal to the stress field in 
the mediimi behind the front. In polar coordinates this is 



a 

rr 



°99 



- - »p - »p 28 

= (^) Op sin 29 
= (^) Op + (^) cos 29 



(B-1) 



is the stress in the mediimi in the direction of wave propagation and vo 
is the stress parallel to the wave front, derived from the assumption 

that there is no strain parallel to the wave front, is equal to 



V 



V 

1-v 



(B-2) 



It can be seen from equations (B-1) that the n = 0 and 2 modes describe 
exactly the free field stresses. Therefore, the unknown stresses at the edge 
of the cavity can likewise be expressed in terms of these modes 



a = - p - p cos 20 
rr o 2 



(b-5) 



a Q = S_ sin 29 
rO 2 

where P^^ P^^ are the unknown modal components of stress acting at 

the boundary. These same stresses must act on the shell. Thus the boundary 
stresses in the medium and shell can be illustrated as shown in Fig. 9* 
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The equations of equilibrium in the medium under plane strain 
conditions are satisfied by the following expressions for the components of 
stress 



_ 1 

"rr = r ar %2 ^^2 

_ ^ ^ 1 

~ 2 dd r ^r^0 

r 



(B-4) 



where Cl is the stress function in terms of r and 0. 



B.ll n = 0 Mode 

The general solution of the stress function for the n = 0 mode is 
given by Timoshenko as 

2 

n = K log r + Mr 

from which can be derived the stress components 



a 

rr 



= ^ + 2M 
r 



a 



r0 



= 0 



(B-5) 



a 



dd 



^ + 2M 
r 



The coefficients K and M are determined from the states of stress at r — > 
and r = R 



K = R 



(¥> % - 



M = - (t^) » 



(B-6) 



For the conditions of plane strain^ the strains are 
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where 





V a. 



rr 



t9 



p °r0 



(B-7) 




modulus of plane strain 
for the medium 



Displacement components are found by suitable integration of the following 
equations for the strains in the medium 



bu. 

^r ^ br 

u 1 

^0 “ r r ^0 

_ 1 ^ ^ V 

'r0 ^ T ho hr r 



(B-8) 



Displacements in the medium for the zero mode are 



u 

R 



E 



a 

P 




V 

R 



0 



(B-9) 



The corresponding displacement components of the shell with an exterior 
compressive force are found in Flugge (2) to be 



u 

s 



R 



P R 
o 

E A 
s 



V 

s 



R 



0 



E 

s 




(B-IO) 



modulus of plane strain 
for the shell 



where 
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The unknown stress is now determined hy equating the radial displacement 
components^ to get 



P 

o 



(l+v) o 



ER,- 

— • X + 1 + V 

.E A 
s 



(B-11) 



After P^ is determined^ stresses in the medium at any radius can he found 
using equations (B-5)« Equation (B«10) gives the displacement of the shell; 
and the hoop and lending stresses of the shell are^ respectively 



P R 
o 

- -j— 

(B-12) 

*^sB Ust 

"FT " 2 R 
' P' 



B.12 n = 2 Mode 



The general solution in terms of a stress function 0, is given as 



Q = 



2 M 

kr + ^ + N 
r 



cos 20 



from which the stress components in the medium become 



rr 



QQ 



t 9 



r r -■ 



cos 20 



PK + ^ 

2K + -^ 
r 



cos 20 



6m ^ 

2K - - -2 

r r ■ 



sin 20 



(B-13) 



Coefficients are determined as Before from the states of stress at r 
and r = R 



K = (ip^) a 

^ 4 p 



(B-lk) 
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M = 



R 



4 r 



T 

N = ^ 
2 



a - P - 2S 

2 p 2 2, 






The displacement components of the medium at the boundary can be expressed 



as 



u 

R 



V 

R 



_1_ 

E 

E 



■ - fS-v) 2(2 -v)S ■ 

2(v-l) + + j— ^ 



2(1-7) . h|r£) 



cos 20 



sin 20 



(B-15) 



and the corresponding components for the shell are 



u 

s 


r ^2 


R^ 


r51 


R 


- E 
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E 



cos 20 



R 



.3 



(B-l6) 



R 



^ , B' R N 2 / R'^ ^ n s 



-E 



sin 20 



The unknown boundary stresses and are determined by equating the dis- 
placement of the shell and medium at the boundary. The hoop and bending 
stresses in the shell are then found from the following equations 



00 



SB 



P R 2S R 
tR3 

(2Pg + Sg) cos 20 



(B-1?) 



Displacements of the shell and stresses in the medium can be determined usxng 
equations given earlier in the discussion. Total static stresses and dis- 
placements are merely the sum of the modal values for any angle. 
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B.2 Shear Wave 

The free field shear wave stresses expressed in polar coordinates 
include only the n = 2 mode 



a = o sin 20 



rr s 

cos 20 (B-18) 

Oaa = -a sin 20 
yy s 

where is the amplitude of the incoming wave. The static solution is 
obtained exactly as in the case of the dilatational wave except for the 
interchange of sines and cosines resulting from the difference in geometry. 
The resulting equations are 

a. Displacement components 



u 

s 

R 



+ 



•9 E„ I l8 E I. 

^ S g 



sin 2G 



V 



'!i- (s! H- H- S)' 

. 6 e ^ V A 5E^ ^121 A . 



cos 20 



h. Shell Stresses 



00 



•P^R 2SgRl 



sin 20 



sB 



IIj ( 2 P 2 + Sg) sin 20 



(B-19) 



(B-20) 



c. Medium Stresses 

a = - (2A + ^ + -^) sin 20 
r r 

= - (2A - ^ cos 20 

r r 

Oqq = (2A + sin 20 

r 



(B-21) 
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where 



A = 



a 

s 



2 



C 




+ ^ ^^2 



D = 





In the machine program^ the n = 0 mode is purely rotational and 
thus we get a static value of tangential displacement. This is calculated 
as follows 



a = 1 y ^ 1 !s. 
rot 2 xy 2 jj, 

o 

1 = 9 = 

R rot 2 [I 



(B-22) 
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DIIATATIOHAL WAVE 



time = 20 R/ci 



Mode 


Quant* 


Stat. 


Mach. 


Dev. Proa 
Static 


Mode 


Quant 


Stat. 


Mach. 


Dev. Proa 
Static 


0 


-oL 


4.000 


4.006 


. 15 ^t 


2 


- OC 


4.613 


4.599 


.30^ 


0 


-<r.e 


4.000 


4.006 


.15 


2 


/3 


3.845 


3.856 


.29 


0 


<Tsb 


.100 


.100 


.00 


2 


Gse 


3.077 


3.114 


1.20 












2 


<Tsb 


. 3^6 


• 3^5 


.29 



* Stresses given in units of |Cp| 

Displaceaents defined by Bqs. (2-50) and (2-36) 



ti»e = 20R/c2 



Mode 


Quant. 


stat. 


Mach. 


Dev. Pros» 
Static 


0 


(3 


5.000 


*^.996 


CO 

0 

• 


2 


OL 


13.840 


13.930 


.65 


2 


/3 


11.535 


U.539 


.03 


2 


“^99 


9.232 


9.149 


.90 


2 


CTsb 


1.038 


1.045 


.67 



XABLE 1 



COMPARISON OP STATIC AND MACHINE VALUES POR MODAL RESPONSE OP THE SHELL 



DILATATIOHAL WAVE 



time = 20 R/ci 







1.0 


R 


0 

01 

i 


R 


3.0 


R 


4.0 


R 


Mode 


Quant? 


Stat. 


Mach. 


Stat. 


Mach. 


Stat. 


Mach. 


stat. 


Mach. 


0 


- Oyy 


.200 


.200 


.550 


.551 


.615 


.616 


.638 


.639 


0 


- Ose 


1.133 


1.135 


.783 


.785 


.718 


.720 


.696 


.697 


2 


- 


-.153 


-.155 


.111 


.112 


.226 


.229 


.271 


.273 


2 


OV® 


.307 


.312 


.i^28 


.429 


.384 


.386 


.363 


.366 


2 


OgB 


.872 


.868 


.367 


.368 


.340 


.342 


.335 


.337 



* Stresses given in units of lffp| 



SHEAR WAVE tisie = 20 R/c2 







1.0 R 


2.0 R 


3.0 R 


4.0 R 


Mode 


Quant. 


Stat. Mach. 


Stat. Mach. 


Stat. Mach. 


Stat. Mach. 


2 


Oyir 


-.460 -.455 


.332 .331 


.675 .678 


.814 .815 


2 




.921 .913 


1.283 1.286 


1.151 1.154 


1.090 1.092 


2 


-Oee 


2.616 2.636 


1.101 1.104 


1.020 1.022 


1.006 1.009 



TABLE 2 



COItPARISOM OF STATIC AlID HACHIHE VALUES FOR MODAL STRESSES HI THE MEDILM 



80 



\ 




(90“) 


DIZ* 


tnf 




(90- ) 


DLF 


t« 


Static 


Max. 


Static 


Max. 


.01 


-9.21 


-10.12 


1.10 


3.12 


-2.49 


-2.73 


1.09 


3.13 


.05 


-7.08 


- 7.74 


1.09 


2.88 


-2.01 


-2.19 


1.09 


2.88 


.10 


-5.56 


- 6.15 


1.10 


2.25 


-1.67 


-1.84 


1.11 


2.25 


.15 


-4.60 


- 5.21 


1.13 


2.25 


-1.45 


-1.65 


1.13 


2.25 


.20 


-3.94 


- 4.52 


1.14 


2.25 


-1.30 


-1.50 


1.15 


2.25 



\ 


<3»b (90“) 


DIF 


tm 




(90”) - 


DIF 


ta 


Static 


Max. 


static 


Max 


.01 


.10 


.11 


1.11 


3.38 


- .09 


- .10 


1.10 


3.25 


.05 


.45 


.50 


1.11 


3.25 


- .35 


- .39 


1.11 


2.88 


.10 


.83 


.94 


1.13 


3.00 


- .55 


- .63 


1.14 


2.50 


.15 


1.18 


1.36 


1.15 


2.88 


- .68 


- .79 


1.17 


2.50 


.20 


1.50 


1.76 


1.18 


2.88 


1 

• 

\Ji 


- .91 


1.21 


2.50 



\ 


oC 


(0“) 


DIF 


tm 


/3 


(45-) 


DIF 


tffi 


Static 


Max. 


Static 


Max. 


.01 


-9.65 


-10.64 


1.10 


3.12 


4.67 


5.23 


1.12 


3.75 


.05 


-8.61 


- 9.50 


1.10 


3.00 


3.85 


4.33 


1.13 


3.50 


.10 


-7.76 


- 8.62 


1.11 


2.76 


3.33 


3.82 


1.15 


3.25 


.15 


-7.15 


- 8.01 


1.12 


2.88 


3.02 


3.55 


1.18 


3.25 


.20 


-6.66 


- 7.56 


1.13 


2.88 


2.80 


3.37 


1.20 


3.38 



* Ratio of Max. to Static # Time in transit times at vhlch Max. occurs 

TABLE 3 



EFFECT OF RELATIVE IHICKHESS OF SECT.! 



31 



\ 


098s 


( 90 ") 


DI2^ 




0e9>« 


( 90 ") 


1 

DLF 


1 

tffi 


Static 


Max. 


static 


Max. 


2 


- 4.16 


- 4.57 


1.10 


2.76 


-2.27 


-2.51 


1.11 


2.76 


k 


- 7.08 


- 7.74 


1.09 


2.88 


- 2.01 


-2.19 


1.09 


2.88 


6 


- 9.33 


- 10.11 


1.08 


2.76 


- 1.81 


-1.97 


1.08 


2.76 


10 


- 12.57 


- 13.41 


1.07 


2.76 


-1.55 


-1.65 


1.07 


2.76 


16 


- 15.72 


-16.47 


1.05 


2.88 


-1.31 


-1.37 


1.05 


2.76 



\ 





DIF 


to 




( 90 “ ) 


DIF 


to 


Static 


Maa. 


Static 


Max. 


2 


.23 


.26 


1.12 


3.25 


- .21 


- .23 


1.13 


3.12 


4 


.45 


.50 


1.11 


3.25 


- .35 


- .39 


1.11 


2.88 


6 


.64 


.71 


1.11 


3.12 


- M 


- .51 


1.10 


2.76 


10 


1.01 


1.12 


1.10 


3.12 


- .63 


- .68 


1.08 


2.62 


16 


1.53 


1.68 


1.09 


3.12 


- .78 


- .83 


1.06 


2.62 



\ 


ct 


( 0 ") 


DIF 


t. 


a 


(»» 5 *) 


DIF 


to 


Static 


1 Max. 

J 


Static 


Max. 


2 


- 4.60 


- 5.12 


1.11 


3.00 


2.14 


2.43 


1.13 


3.62 


4 


- 8.61 


- 9.50 


1.10 


3.00 


3.85 


*^.33 


1.13 


3.50 


6 


- 12.22 


-13.37 


1.09 


2.88 


5.32 


5.96 


X .12 


3.38 


10 


-18.65 


- 20.20 


1.08 


2.88 


7.92 


8.82 


1.11 


3.25 


16 


-27.11 


-29.15 


1.08 


3.00 


11.41 


12.64 


1.11 


3.12 



*Batlo of Max. to Static In transit tlaes at which Meuc. occurs 

TABLE k 






EFFECT <»■ RATIO OF MODULI OF ELASTICm 



82 



> 




(90") 


DLP* 






(90*_) 


DU 


tm 


Static 


Max. 


static 


Max. 


1 


-7.08 


-7.67 


1.08 


3.00 


-2.01 


-2.17 


1.08 


2.88 


3 


If 


-7.7^ 


1.09 


2.88 


If 


-2.19 


1.09 


2.88 


5 


It 


-7.97 


1.13 


2.50 


tl 


-2.26 


1.12 


2.50 


7 


tl 


- 8.30 


1.17 


2.61 


It 


-2.3*^ 


1.17 


2.62 


10 


11 


-8.70 


1.23 


2.88 


If 


-2.45 


1.22 


2.88 



> 


Osb 


(90’I. 


DU 


tja 




(90") 


DU 


1 

to 


static 


Max. 


Static 


Max. 


1 


.45 


.49 


1.10 


3.25 


- .35 


- .38 


1.09 


3.00 


3 


H 


.50 


1.11 


3.25 


H 


- .39 


1.11 


2.88 


5 


It 


.50 


1.13 


3.12 


It 


- .41 


1.15 


2.88 


7 


tl 


.51 


1.15 


3.25 


It 


- .^3 


1.21 


3.12 


10 


It 


.53 


1.18 


3.50 


tl 


- .46 


1.30 


3.50 



\ 


OL 


(0*1. 


DU 




f3 


(45") 


DU 


tm 


Static 


Max. 


Static 


Max. 


1 


-8.61 


-9.38 


1.09 


3.00 


3.85 


4.27 


1.11 


3.50 


3 


tf 


-9.50 


1.10 


3.00 


M 


4.33 


1.13 


3.50 


5 


H 


-9.64 


1.12 


2.88 


tl 


4.4l 


1.15 


3.62 


7 


H 


-9.77 


1.14 


3.00 


11 


4.45 


1.16 


2.76 


10 


« 


-9.95 


1.16 


2.62 


It 


4.72 


1.23 


3.00 



*Ratio of Max. to Static ^^^Tiae in transit tiaes at which Max. occurs 

TABLE 5 



EEPBCT MASS RATIO 



83 



t) 




( 90 *) _ 


DU** 


tj 




( 90 °) 


DLF 


tm 


Static 


Max. 


static 


Max. 


0 


- 8.33 


- 9.17 


1.10 


3.25 


-2.08 


-2.29 


1.10 


3.25 


.1 


- 7.96 


- 8.69 


1.09 


3.12 


-2.06 


- 2.24 


1.09 


3.12 


.2 


- 7.42 


- 8.08 


1.09 


3.00 


- 2.02 


- 2.20 


1.09 


3.00 


.25 


- 7.08 


- 7.74 


1.09 


2.88 


- 2.01 


-2.19 


1.09 


2.88 


.4 


- 5.76 


- 6.75 


1.17 


3.12 


-1.91 


- 2.22 


1.16 


3.25 



V 


<Tsb 





DI^ 


tm 


Ovv yy\ 


( 90 °) 


DLF 


tn 


static 


Max. 


static 


Max. 


0 


.61 


.69 


1.14 


3.25 


- .42 


- .47 


1.13 


3.00 


.1 


.56 


.63 


1.13 


3.25 


- .40 


- .44 


1.12 


2.88 


.2 


.49 


.55 


1.12 


3.25 


- .37 


- .hi 


1.11 


2.88 


.25 


.45 


.50 


1.11 


3.25 


- .35 


- .39 


1.11 


2.88 


.4 


.27 


.31 


1.14 


3.75 


- .29 


- .33 


1.15 


3.88 





(X 


( 0*1 _ 


Dir 


tn 


<3 


(45* )_ 


mr 


tn 


static 


Max. 


static 


Max. 


0 


-10.33 


- 11.48 


1.11 


3.25 


6.00 


6.94 


1.16 


3.38 


.1 


- 9.88 


-10.90 


1.10 


3.12 


5.31 


6.08 


1.14 


3.38 


.2 


- 9*12 


-10.03 


1.10 


3.00 


4 . 4 o 


4.98 


1.13 


3.38 


.25 


- 8.61 


- 9.50 


1.10 


3.00 


3.85 


4.33 


1.13 


3.50 


.4 


- 6.54 


- 7.61 


1.16 


3.50 


1.77 


1.97 


1.11 


4.38 



Ratio of Meoc. to Static "Tlae in transit tines at which Max. occurs 



TABUS 6 



EFFECT OF POISSOH'S RATIO 




Plaae Dilatatlonal Wave Front at time t - 0 
Stresses Behind Wave Front Shown In Both Rectauogular 
and Polar Coordinates 




Plane Dilatatlonal Wave Front at time t > 0 



FIGURE 1 




Plane Shear Wave Front at time t = 0 
Stresses Behind Wave Front Shown in Both Rectangular 
and Polar Coordinates 




Plane Shear Wave Front at time t > 0 



FIGURE 2 
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FIG. 3 DISPLACEMENT COMPONENTS OF THE SHELL 




FIG. 4 DEFINITION OF POSITION ANGLE, G| 
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FIG. 5 GZCKSm or A POIRT SOURCE ON OBB z AXIS 




PIG. 6 GBOMETRI OF MANY POINT SOURCES 




PIG. 7 POSITION ANC0LE, 0, , FOR ANY RADIUS 




FIG. 8 TIME DEPENDENT STRESS WAVE 





FIG. 9 SSATIC STRESSES ACTING ON IBE MEDIUM AND SHEIX 
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FIO. 10 WEKHtTIHO COEFFICIEHTS Al^ and BMn, 
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PIG. 11 FUNCTIONS PC-I^T) and 0(-l/lCc+‘«' ) 
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FIO. 12 FUNCTIONS and 0(-l/kc^‘r) 
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FIG. 13 COMPOHENTS OF ACCELERATICaf FOR MODE n * 0 and 1 
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FIG. 15 MODAL VELOCITI C< 



96 




MD l/N ^ or> OJ rH O 



00 



VO 



ITN 




fO 



CVi 



O 



\*io\ jD SLiMn 



FIG. l6 MODAL DISPLACEMENT CCMPONENTS 
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FIG. 17 MODAL CCMPOKENTS OF SHELL STRESSES 
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TRANSIT TIMES 

FIG. 18 MODAL COMPONENTS OF MEDIUM STRESSES 
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PIG. 19 MODAL DISPlACEMEarr COMPONENTS FOR INCUEHT SHEAR WAVE 
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Fia. 20 MODAL SHELL STRESSES FOR INCIDENT SHEAR WAVE 
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PIG. 21 HOOP AND BEMDIHO STRESSES THE SHELL AT VARIOUS AKOI£S 
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FIG. 22 RADIAL AMD SHEAR STRESSES IM THE MEDIUM AT THE 



103 




CO 






VO 



l/^ 




on 



OJ 



rH 



o 



0 1^0 O O 

« • • • • • • 

CO CM CM fH H o O 



I*!£>I io sxDin 



PIO. 23 HOOP STRESSES IH THE MBOHM AT THE 
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FIG. 2k STBS8SS3 IN THE SHELL AND MEDIIM FCR INCIDEINT SHEAR WAVE 
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FIG. 25 HOOP STRESSES IN THE SHELL AND MEDIUM FOR VARIOUS RELATIVE THICKNESSES OF SHELL 
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FIO. 26 RADIAL STRESS IN THE MEDIUM AT THE BOUNDARI FOR VARIOUS THICKNESSES OF SHELL 



107 




FIO, 27 SHEAR STRESS IN THE MEDIUM AT THE BOUNDARY FOR VARIOUS RELATIVE THICaOIESSES CF SHELL 
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FIG. 28 STRESSES IN IHE MEDIUM FOR VARIOUS RADII AT TEN TRANSIT TIME 





PIG. 29 VARIATION OF STRESSES IN THE MEDIUM WITH RADIUS AND TIME 
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FIG. 30 STRESSES IN THE MEDIUM FOR VARIOUS RADII AND TIMES 
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TRANSIT TIMES 

FIG. 31 CCMPARISON OF SHORT-TIMB RESULTS FOR INCIDENT P WAVE 
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FIG. 32 COMPARISON OF SHORT-TIME RESULTS FOR INCIBENT P WAVE 
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TRANSIT TIMES 

FIG. 33 COMPARISON OF SHORT-TIME RESULTS FOR INCIIENT P WAVE 
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TRANSIT TIMES 

PIG. 34 CCMPARISON OF SHORT-TIME RESULTS FOR INCIDEHT S WAVE 
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FIG. 35 COMPARISON OP SHORT-TIME RESULTS FOR INCIIffiNT S WAVE 
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FIG. 37 EFFECT OF SIZE OF TIME INTERVAL ON LONG-TIME RESULTS 
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PIG. 38 EFFECT OF RELATIVE THICKNESS OF SHELL 
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FIG. 39 EFFECT OF RATIO OF MODULI OP ELASTICITT 
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PIG. 40 EFFECT OF MASS RATIO 
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FIG, 4l EFFECT OP POISSC»J'S RATIO 
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FIG. 42 EFFECT OP ADDITICaiAL MASS (W THE RESPONSE OF THE SHELL 



123 




00 



\o 



ur^ 




m 



CVI 



o 



VO 






1^1 JO siiun 



CVI 



o 

H 



00 



CM 



O 



PIG. U3 EFFECT OF TIME VARYING INCIDENT STRESS WAVE ON SHELL HOOP STRESS 
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FIG. 44 EFFECT OF TIME VABYIHO IWCIDENT STRESS WAVE GS SHELL HOOP STRESS 
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PIG. 45 COMPARISON WIIB PREVIOUS WORK, HOOP STRESS (SHELL) FOR n 
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FIG. U6 COMPARISON WITH PREVIOUS WORK, HOOP STRESS (SHELL) FOR 
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PIG. 47 COMPARISON WITH PREVIOUS WORK, RADIAL DISPLACEMENT FOR n 
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FIG. 48 COMPARISON WITH PREVIOUS WORK, RADIAL DISPLACEMENT FOR n 
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FIO, h9 COMPARISON WITH PREVIOUS WORK, TANGENTIAL DISPLACEMEWT FOR n 
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